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Abstract. We reconsider the functional renormalization-group (FRG) approach to decaying Burgers turbulence, 
and extend it to decaying Navier-Stokes and Surfacc-Quasi-Geostrophic turbulence. The method is based on a 
renormalized small-time expansion, equivalent to a loop expansion, and naturally produces a dissipative anomaly 
and a cascade after a finite time. We explicitly calculate and analyze the one-loop FRG equations in the zero- 
viscosity limit as a function of the dimension. For Burgers they reproduce the FRG equation obtained in the 
context of random manifolds, extending previous results of one of us. Breakdown of energy conservation due to 
shocks and the appearance of a direct energy cascade corresponds to failure of dimensional reduction in the context 
of disordered systems. For Navier-Stokes in three dimensions, the velocity-velocity correlation function acquires a 
linear dependence on the distance, £2 = 1, in the inertial range, instead of Kolmogorov's £2 = 2/3; however the 
possibility remains for corrections at two- or higher-loop order. In two dimensions, we obtain a numerical solution 
which conserves energy and exhibits an inverse cascade, with explicit analytical results both for large and small 
distances, in agreement with the scaling proposed by Batchclor. In large dimensions, the one-loop FRG equation for 
Navier-Stokes converges to that of Burgers. 
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1. Introduction 

Describing Navier Stokes (NS) turbulence with the tools of statistical physics has remained a major challenge 
since Kolmogorov's dimensional arguments leading to the E{k) ~ ^2/3^,-5/3 cncr gy spectrum for the 3D energy 
cascade [1, 2, 3, 4]. The simplest analytical method, Kraichnan's direct interaction approximation closure scheme [5] 
(equivalent to mode coupling) failed to recover Kolmogorov's prediction. There were numerous attempts to overcome 
these difficulties using a variety of methods, e.g. more refined closure schemes [6], large number of components [7, 8], 
renormalization-group (RG) [9, fO, 11, 12, 13, 14, 15] conjectures for short-distance expansions [16, 17, 18], study of 
short time singularities [19], tetrad models [20], and shell models [21], with various degrees of success. At the heart 
of the cascade phenomenon is that non smooth velocity field do not conserve energy. The main physics challenge, 
i.e. to describe the statistics of the energy transfer via singular or almost singular structures, is only partially 
captured. Predicting the multi-fractal corrections for velocity moments S p (u, t) = ([(v ut — vot) • u/|u|] p ) ~ C p |u|^ 
to Kolmogorov's prediction (C^ = p/3) remains a challenge, despite the analytical progress in the simpler passive 
scalar problem [22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40]. For the inverse cascade 
in 2D [41], due to an infinity of conserved quantities, and a simpler numerical modeling, more is known; the most 
recent analysis unveils a tempting connection to conformal field theory and SLE [42] , but remains based on numerics 
or speculative [43]. Also relations between Navier-Stokes and the membrane dynamics of black holes have been 
discussed in the context of the AdS/CFT correspondence [44, 45, 46] as well as relations to the physics of graphene 
[47]. 

The problem of ./V-dimensional Burgers turbulence, i.e. of a potential flow without pressure, exhibits similarities 
with Navier-Stokes turbulence, such as the existence of an inertial range which supports an energy cascade and the 
multi-scaling of the velocity moments. Although, as NS, it lacks a small control parameter and hence is non-trivial, 
it is simpler, since the Burgers equation can be integrated explicitly via the Cole-Hopf transformation J, hence it has 
allowed for some progress [50, 48]. A remarkable mapping to an elastic object in a quenched random potential maps 
the shocks in both a decaying or stirred Burgers velocity field to the jumps of the equilibrium position of the pinned 
elastic object (which is a point for decaying Burgers or a line for stirred Burgers) upon variation of an external field. 
This mapping was used to study the large-dimension N limit of stirred Burgers turbulence using replica symmetry 
breaking [51, 52] and, more recently, of decaying Burgers turbulence [53, 54]. The detailed statistics of shock cells 
which is obtained from these works is consistent with the physical expectation, and important open questions are 
now (i) whether this is a good starting point to perform an expansion towards finite TV; (ii) whether it can inspire an 
approach to Navier Stokes in large dimension, a notably difficult problem [55]. An RG inspired method bypassing 
the Cole-Hopf transformation has been proposed very recently for the KPZ equation which is closely related to the 
Burgers equation [56]. 

Another powerful method able to handle singularities such as shocks and avalanches in disordered systems, 
which does not rely on large N , is the Functional Renormalization Group (FRG) [57, 58] (for an introduction and 
review see [59, 60]). The connection between the FRG and decaying Burgers turbulence was elucidated in [61, 62] 
(see [63] for an earlier attempt). It turns out that the force felt by an elastic manifold of internal dimension d 
submitted to a random potential plus a quadratic well can be seen as a generalized velocity field: it satisfies an 
exact evolution equation which is a functional generalization of the decaying Burgers equation, where the role of 
time is played by the (inverse) curvature of the well. For d = the manifold is a point and one recovers the standard 
Cole-Hopf representation of the Burgers equation. The hierarchy of equations relating 71-point equal-time velocity 
correlation functions identifies with the (exact) hierarchy of FRG flow equations, and the loop expansion in the 
field theory corresponds to the (rcnormalizcd) small time expansion in the (generalized) Burgers problem, as will be 
detailed below. The amazing property is that this hierarchy becomes controlled in an expansion in e = 4 — d around 
d = 4, which is the crucial property of the FRG approach to disordered systems. Hence Burgers turbulence, i.e. 
d = 0, becomes accessible via this expansion. Furthermore the physics of the generalized Burgers problem (i.e. of 
the manifold) has features which are independent of the parameter d. For instance, energy conservation for smooth 
flows is obtained as well as an infinite number of conserved quantities (the first property being called " dimension 
reduction" in the context of manifold, and the second corresponds to the non-rcnormalization of the moments of 
the so-called Larkin random force). Non-conservation of energy via shocks occurs for any d, and the dissipative 
anomaly at the heart of the energy cascade, i.e. the non-vanishing limit of the energy flux —dt£ = e = v((Vv) 2 ) 
as v — > 0, is naturally captured by the FRG [61, 62]. Finally the FRG allows to compute shock-size distributions 
in the controlled expansion around d = 4 [64, 65]. Most of these studies focused on N = 1, but recently we also 

f It is often said [48] that it lacks an essential property of Navier-Stokes turbulence, namely the sensitivity to small perturbations in the 
initial data, and thus the spontaneous appearance of randomness by the chaotic dynamics. However, there is, in some cases (denoted 
SR below), decorrelation in time of two slightly different initial conditions, a property sometimes termed chaos in the community of 
disordered systems [49]. 
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investigated N > 1 [66]. 

The aim of this paper is thus to investigate whether FRG-inspired methods can be developed to describe 
Navier-Stokes turbulence as well. Here, our scope is relatively modest and it should be seen as a first exploration 
of the FRG method into the domain of non-linear physics. We focus on the decaying Burgers and Navier-Stokes 
equations; however, the stirred case can also be studied within the same framework. We derive the 1-loop FRG 
equations, first for Burgers in A-dimcnsion (since most explicit calculations in [61, 62] were for N = 1), and then 
for Navier-Stokes. We discuss some features of the fixed-point solutions which correspond to a decaying turbulent 
state, leaving a detailed analysis for the future. In N — 2 we perform a detailed analysis of the NS fixed point. 
At this stage, the method for NS is not a controlled perturbative expansion scheme, since there is no equivalent of 
the Cole-Hopf mapping. The method however does capture some of the physics of the singularities. We analyze 
the nature of the singularities at small distance. While our analysis is restricted to one loop, we discuss at the end 
possible extensions to higher loops. 

Let us stress that in the so-called inertial range of length and time scales, it is widely expected that the statistics 
of decaying turbulence is rather similar to the forced one. Indeed, due to strong separation of the large and small 
time scales, the eddies in the inertial range have enough time to reach an equilibrium for the energy-flux, before 
the larger eddies will significantly decay The scaling behavior of decaying turbulence in the inertial range is thus 
indistinguishable from the forced case, while differences will occur at the scales of the large eddies. This is the 
universality assumption entering most theories of turbulence, see e.g. Ref. [67] for a detailed discussion. 

The paper is organized as follows: In section 2, we introduce the model and notations. In section 3, we review 
known results, both for Burgers and Navier Stokes. The FRG equations are derived in section 4. We start by the 
general scheme, before giving results for Burgers and Navier Stokes; finally we discuss conserved quantities. In 
section 5, we discuss the short-distance singularity of the FRG equations: Are there solutions with other power 
laws than a cusp? In section 6 we focus on the analysis of the FRG equations for Navier-Stokes turbulence in two 
dimensions. In section 7 we discuss the limit of large N and in section 8 we give the FRG equations for surface 
quasi geostrophic turbulence. 

Finally note that this work was started a while ago. For an early exposition sec [68]. 



2. Model and notations 



We study here two models: 

(i) the A-dimensional decaying Burgers equation for a A-componcnt velocity field v ut at point u and time t, 

d t xZt = vVKt - \dZ{^tf • (1) 

The velocity is assumed to be vorticity-free so it can be expressed as gradient of a potential function, = <9„ V(u, t), 
that implies ^d a (v) 2 = (v ■ V)v Q . Note that to streamline notations we attach space and time indices to the fields, 
which are never to be understood as derivatives. We use boldface to indicate vectors, and normal font for scalars, 
so k := |k|. 

(ii) the incompressible A-dimensional Navier-Stokes equation 

e t vz t = - pjp(d u )d2(viy ut ), p^(d u ) = - (2) 

with the pressure eliminated using the transverse projection operator P^g(d u )- The latter implies the divergence 
free constraint (incompressibility) V • v = at all times. In Fourier space, both equations can be written as 

dtvZ t = -vk\Zt ~ ^*Wk) £ yfrk, (3) 

p+q=k 

with 

P 3 fk^-J ikaS ^ (Burgers), 
a '^ [ ' ~ \ tk*P£(k) + ikyP^ (k) (Navicr Stokes). 1 > 

The transverse and longitudinal projection operators written in Fourier space read 

In both cases we are interested in the small-viscosity (large Reynolds number) limit v —> 0, in which case a broad 
inertial range develops. In that limit the Navier-Stokes equation formally becomes the Euler equation, and in both 
cases weak solutions exist [69, 70, 71] (for review see [72]). 



W = ^-^> p Q yk ) = ^. ( 5 ) 
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We study homogeneous turbulence with random initial conditions, i.e. we chose an initial condition at time t = 
which is statistically translational invariant. Everywhere we denote (. . .) the average over initial conditions. We 
assume that the initial velocity field is Gaussian distributed and that its spectral support is around a characteristic 
wave number fco. The averaged squared initial velocity is vq = (vq ) 1//2 . The initial Reynolds number is 
1Z = Vo/(uko) which wc assume to satisfy 1Z ^> 1. The initial range (where viscosity is subdominant), given 
by kg <C k <C vq/v, has no particular structure, but will develop, as time increases, into a self-similar inertial range 
where the energy cascade (in three dimensions) will take place. 

To stress the similarity with the FRG, we denote the equal-time velocity two-point function as 



'ut v 0t 



g t ) = A w (u) , <v£ t v£, t ) = a k - k ,A t , Q/J (k) . (6) 

We denote := (2tt) n S N (k), and in Fourier space A{ >a| g(k) = P a i (k)At(k) for Burgers and Navier-Stokes 
respectively. By definition, A Q( g(u) = Ap a (— u), while symmetry is only assured for isotropic turbulence. We study 
a system in a periodic cube (torus) of volume L N with mostly two distinct cases: 

(A) fixed t and L — > oo in which case we further restrict to isotropic (homogeneous) turbulence where A t (k) 
depends only on |k|. 

(B) fixed L, which becomes periodic turbulence. In case (B) we use discrete Fourier modes 5Z q , which implicitly 

become J2q Jq = / jhrW m an formula below in case (A). 

The total kinetic energy per unit volume, and the kinetic-energy spectrum arc denoted as in [6] 

1 r°° 

£(t) = 2 {v ' 2)= J dfcf (M) • (7) 

However for convenience we also use a non-standard normalization for energy and energy spectrum and denote (Sn 
is the area of the unit sphere in dimension N): 

E(t):= [ —^(v 2 )= / dkE(k,t) , (8) 

Jo 

E(k, t) := J —5(k - |k|)A 4 , Qa (k) , (9) 

which for isotropic turbulence becomes 

E(k,t) =k N - 1 A t , aa (k) . (10) 

Hence £(k,t) = \ ( 2 j"« E{k, t). The decay of turbulence depends on the initial spatial behavior of the energy 
spectrum at large scales, hence we denote (in the isotropic case) 

E(k,t = 0) ~ fe _ ck n , A , QQ (k) ~ k ^ Q ck n - N+1 , A , QQ (u) ~ U _ KX) c'w-"- 1 . (11) 

The total (kinetic) energy (per unit volume) E(t) = dkE(k 7 t) ~ A tjQ , a (u = 0) must be finite (it grows as L n+1 
for n < —1) and the initial velocity field is usually assumed to be smooth, i.e. Ao, a ^(u) is an analytic function in 
each of its components u near |u| = 0. 

In all cases below, when the system reaches a statistically scale invariant decaying state, we denote by 

t(t) = t C/2 (12) 

the characteristic length scale, which usually separates the inertial range from the large-scale region (i.e. i(t = 0) ~ 
1/ko). The notation is motivated by the relation, in the Burgers case, to the roughness exponent £ for random 
manifolds. We distinguish it from the exponent (2, which describes the leading short-distance singular behavior of 
the two-point equal-time velocity correlator in the inviscid limit. We will also assume a dissipation scale £d(t) <C £(t) 
which is the lower boundary of the inertial range, and is often set to zero in the following, corresponding to the 
inviscid limit. 



3. Known results and phenomenology 

3.1. Decaying Burgers 

In the following we use the mapping of the Burgers equation onto the problem of a particle in the iV-dimensional 
potential W = V(uq) + ^j(u — Uo) 2 , where V(u) is the random potential which parameterizes the initial condition, 
i.e. v u t = o = V u V^(u). Denoting by Z the canonical partition function of the particle at temperature T = 2v, the 
velocity at all times is v ut = V u V r (u, t), where V^u, t) = —Tin Z is the free energy. In the inviscid limit it becomes 
V(u,t) = min uo W. Let us summarize what is expected for Burgers (most results are shown for N = 1, and 
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conjectured for N > 1), for more details see e.g. [48] and the discussion in Section E of [62] including connections 
to the FRG. We assume a smooth Gaussian initial velocity field with the correlator (11), i.e., with a spectrum 
proportional to \k\ n at small k and decreasing quickly at large k. For isotropic turbulence (A) there are two cases: 

(1) Long-range initial condition (LR) n < 1: The correlator of the random potential ([V(u) — V^O)] 2 ) grows 
as at large u and the particle is always in a glass phase, i.e. the effective viscosity scales to zero (see below). 
The evolution is expected to reach an asymptotic statistically scale invariant form v ut = £(t)t v(u/£(t)) (in 
law) with £(t) = < c / 2 , where C = 4/(3 + n) and with energy decay E(t) ~ t" 2+c = t-2(«+i)/(n+3)_ Shocks, i.e. 
codimension-one manifolds (together with some additional lower dimensional singularities for N > 1), where the 
velocity is discontinuous (at v = + ) or nearly discontinuous (at small v > 0), form in finite time and, convected 
by the flow, keep merging when they meet. The growing scale of this coarsening process (quite complicated for 
N > 1) is expected to scale as £(t) ~ t?/ 2 . This is clear at least for TV = 1 §. While the width of an isolated shock 
grows as L^t) ~ id, the width of the surviving shocks grows as L' d (t) ~ zy£ 1- ?/ 2 [61, 62], hence the rescaled width 
L' d (t)/£(t) ~ vt 1 ^^ = is c g scales to zero for n < 1. This corresponds in the FRG to an attractive zero- viscosity (i.e. 
zero-temperature) fixed point (describing a glass phase for the particle): One can write u c g = 

For a random manifold, the glass exponent is 8 = d — 2 + 2£ > 0, with here d = 0. This is a LR fixed point, with 
exponents given by their dimensional values (also called Flory values in the context of elastic manifolds), i.e. they 
are fixed by the initial condition. This property, called the persistence of large eddies in turbulence, means that the 
energy spectrum for k <C l/£{t) retains its original form (11) with an amplitude c independent of time. 

(2) Short-range initial condition (SR) n > n c = 1: This is the Kida regime [48, 73] with an asymptotic 
statistically scale invariant form (see however below) with a scale £{t) = ^/(lni)] 1 / 2 and a decay of the energy 
E{t) ~ l/^lni) 1 / 2 ] (for gaussian initial conditions). The (rescaled) shock width now grows (there is no glass 
phase); hence it exists only for v — > before t — > oo (see [74, 75] for a more refined analysis of the double limit). 

There is an additional crossover region (e.g. for 1 < n < 2 for N = 1) where the persistence of large eddies 
(i.e. of the tail of the FRG function) still holds, but the system flows to the SR (Kida) fixed point: This is known 
as the Gurbatov phenomenon [76], i.e. the velocity statistics is not scale invariant. The resulting energy spectrum 
then consists of three regions: (i) the "outer region", < k < k s (t) ~ t~ 1 '^ 2 ~ n " , where the velocity correlations 
preserve its initial form (11); (ii) the "inner region", k s (t) < k < £(t) ~ t -1 / 2 (£ = 1), with spectrum k 2 , and (iii) 
the shock-dominated region with spectrum k~ 2 for k > £{t)~ . All scales are given up to logarithmic corrections 
and for N = 1. More details can be found in [76]. Within the FRG analysis, this crossover region can be seen as 
crossover from the LR to the SR FRG fixed point ||. In the FRG analysis of random manifolds a similar crossover 
was described in [77]. 

In the marginal case n = 1, a LR fixed point exists where the potential retains logarithmic correlations, with 
a phase transition as a function of v [75] . 

Note that originally Burgers [78] distinguished only two cases (for N = 1), assuming that J = A(k = 0) = 
J duA(u) exists. The case J > then corresponds to the LR case n = 0, hence £{t) ~ t 2 ^ (£ = 4/3) and is usually 
called "random-field" fixed point in the language of random manifolds. The case J = was solved by Kida [73], 
and corresponds here to the SR case n > 1 (for instance for n = 2, A(fc) is analytic and the random potential 
has (^-correlations). This is usually called the "random-bond" fixed point in the language of random manifolds. 
The summary presented above contains many more cases, i.e. the LR models form a line of FPs, continuously 
parameterized by n, and the SR case can also be modified by the Gurbatov LR-SR crossover. 

Finally for the periodic case, the system converges, for N — 1, to a single random shock per period with 
E(t) ~ t~ 2 . This corresponds to the FRG random-periodic FP f = (i.e. n = oo), and a similar picture should 
hold for any N. 

3.2. Decaying Navier Stokes 

A similar discussion can be given for decaying NS, though on a much less firm basis, mostly phenomenology, scaling 
arguments, closure calculations and some support from experiments. Again since Von Karman and Howarth [79] one 
assumes a decaying state v ut = £(t)t~ 1 v(u/£(t)) (in law). Then v satisfies a equation where v — !> vtj£(t) 2 , which 
flows to zero if C > 1- In Fourier this can be written v^i = £(t) 2 t~ 1 -v(k£(t)') (in law) and the energy spectrum takes 
the form E(k, t) = £{tf /t 2 E(k£(t)) with £(t) = t c/2 and a total kinetic energy decay E{t) - t~ 2+c . The persistence 
of large eddies, i.e. the invariance of E(k, t) ~ ck n , implies ( = 4/(3 + n). This corresponds to a LR initial condition 
(i.e. regime (1) in the previous section). The energy spectrum can then be divided into a low-wavenumber range 

§ At least for dilute shocks n < — 1; shocks are dense for n = 2 and the analysis is more difficult. 

|| For instance the value n = 2 corresponds to the Flory exponent £ = 4/5, and at short scale the SR correlator of the random potential 
is behaves effectively as 8(x) ~ 1/x, while at large scale it flows to the SR Kida FP. 
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k£(t) < 1 with E(k,t) ~ cfc", and the inertial range k£(t) > 1, with, in 3D (N = 3), £(k,t) ~ C K e(t) 2 / 3 k- 5 / 3 
(assuming Kolmogorov spectrum) with dE/dt = —e(t). For more details see [6] VII- 10, Ref. [80] for RG arguments, 
and [81]. 

There is agreement that this LR regime cannot hold for n > 4 because of the C(t)k 4 backtransfer in the energy 
spectrum discovered by Proudman and Rcid [82] for the 3D Navier-Stokes equation, and found in EDQNM closure 
calculations [6] (analogous to the k 2 backtransfer for the Burgers dynamics). In other words, the low-fc energy 
spectrum cannot be softer than fc 4 . For n > 4 it is argued that the small- k part of the spectrum is replaced by 
E(k,t) ~ C(t)k A at small k, with the inertial range at large k. Because of the Gurbatov phenomenon (analogous 
to the situation in Burgers discussed above) it is then argued [83] that the LR regime cannot hold for n < n c with 
3 < n c < 4. In the range n c < n < 4 there are three spectral regions E(k,t) ~ ck n for k < l/£*(t) (outer region), 
E(k,t) ~ C(t)k A for l/t t (t) < k < l/£(t) (inner region) and finally the inertial range for k > 1 /€(£), leading to a 
breakdown of global self-similarity. This global picture seems compatible with experiments [84, 85]. 

Note that in 3D NS there is another conserved quantity, the helicity h = e Q ^ 7 v Q ,9 / gv 7 . It is locally fluctuating 
even if its average is zero. If its average is non-zero, as in MHD, then we need to consider A aj g(u) ^ A J g Q (u) = 
A a p(—u). Its presence makes possible a joint cascade with two fluxes, one of energy and one of helicity, both to 
small length scales [86]. 



4. FRG equations 



4-1- Loop expansion: General strategy 

We now write FRG-like equations able to access directly the strong-coupling regime (i.e. finite non-linearity) using 
either a graphical method or, equivalcntly, starting from the exact infinite hierarchy obeyed by the equal-time 



n-point correlation functions denoted here (v";,v" 2 



Ult *U2f 



dtCM 



(ui,u 2 ,...,u„) =nu Sym 



• 2 q(t, 



^Sym Pa li ^(V Ul )C , ££i?.. a 



(ui,u 2 , ...,u 
ln (ui,u 2 , ...,u 

(ui,Ui,U 2 ,...,U„) 



They obey, for n>2, 



(13) 



The time dependence is implicit, and -P Q; /3 7 (V) is given by Eqs. (4) and (5) rewritten in real space using zk 4-> V; 
symmctrization with respect to the n pairs (uj, cti) for i = 1, ...,n is denoted by Sym[...]. We recall that 
Cqiq 2 (u x , u 2 ) = A tjQ1Q2 (ui — u 2 ). In Fourier space the hierarchy reads 



5tCW.a n (k 1 ,k 2 ,..,k n )= -ni/Sym k?cW. aiB (ki,k a ,...,k 



Sym 



(n+l) 



\ v kif" v k„t/ 



(p,q, k 2 ,.. 

Cai . ..(X n 

(k!,...,k n )J 

ki+...-t-k n i 



k n ) 



p+q=ki 

V(n) 



(14) 



and we also define 



for the correlations 
C* QlQ2 (k,-k) = C Ql " 2 (k). 

If equation (13) is considered in the inertial range, i.e. all 3> id, it is expected (and for Burgers in some 
cases shown) that one can neglect the viscosity term in the hierarchy. To study the inertial range it is thus tempting 
to consider the limit v = + of these equations. To argue that this can be done, and that the result is still given 
by equation (13) setting v — 0, we need two conditions: (i) the physical requirement that the n-point velocity 
correlations are continuous functions, i.e. that limits at coinciding arguments exist; (ii) the property that the 
v — > limit of averages such as (^(u i )Pj (d u )d^ l (v2 lt 'v^ lt )) where $ is any product of velocities with ^ u is equal 



to P^ /3 (9 u )92( c f ) (u i )v^ t v^ t ). This appears to be related to the existence of weak solutions of the Eulcr equation, 
which is discussed in [69, 70, 71, 72]. In some cases, e.g. for the inviscid Burgers equation and N = 1, it can be 
justified [61, 62] from the dilute shock picture of [87]. 

Assuming that the v = hierarchy holds, let us describe the strategy of the loop expansion as it was constructed 
in the case of Burgers [61, 62]. It amounts to looking for a solution of the hierarchy in the schematic form 
(complicated convolutions are indicated by *): 

d t A = tA * A + t a A * A * A + ... (15) 

C (3) =tA* A + t 3 A* A* A + ... , C (4) = A* A + t 2 A* A* A + ... (16) 

C* (5) =tA* A* A + ... , 



C (6) = A * A * A + . 

It is illustrated here to two loops, and more generally C*-™-* = J2 q >o t 2q+,iq C^ 
even and odd. We impose that at t = the distribution is Gaussian, hence the functionals Cq 2 ^[A], k > 1, are 



(17) 

A] with e q = 0, 1 for n > 3 respectively 
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given by the Wick decomposition. This allows to compute iteratively all the terms in the beta function (15): e.g. 
to one loop we start from Cq = 3[AA], and use (13) with successively n = 3 and n = 2, first to get Cq, then to get 
the one-loop A 2 term in (15). Thus, the beta function appropriate to the rescaling A = t~ 2 A is obtained directly. 
Counter-terms are produced automatically by successive corrections, due to consistent evaluations of dtA terms at 
each step. Higher- loop calculations will be presented elsewhere [88]. The first corrections to each cumulant, i.e. 
the c[ 2k ^ [A] and CQ 2fe+1 ^ [A] arc the tree approximation. For calculations per se, an equivalent procedure, which 
we also performed to one loop using a graphical method directly on the Burgers and Navier-Stokes equations, is 
to compute A t = A t= o + J2 q t 2q ^t=o as a direct small-time expansion, then compute dt A t and re-express the 
result in terms of A t itself as given in (15) by inverting the series. The information contained in the beta function 
can thus be described as a renormalized small-time expansion (for a direct small-time expansions see [55, 19]). 

Note that the present FRG is different from the usual RG for turbulence as developed in [9, 10, 11, 12, 13, 14, 15]. 
Here we keep the complete crucial information about A(k) (in Fourier) while the usual RG integrates out shells in 
k. In particular, the information in the small-u behavior of A(u) ~ iS 2 determines the exponent £2. 



.2. FRG equation for Burgers 



We now display the resulting one-loop /3-function for Burgers. Calculations in real space are given in Appendix A, 
and in Fourier space in Appendix B. A graphical derivation can be found in Appendix D. The result reads % 

2 



dtA aP (u) = -td a df3 A 75 (u) - A 75 (0) 



(18) 



It has the usual form of the A-component one- loop FRG equation [77, 89] . To study the evolution from the initial 
condition (11), it is more convenient to introduce the rescaled velocity correlation A tiQ( g(u) defined through 



A t , Q(3 (u) = i< c -4At, Q(3 (u/^ /2 ) . 



(19) 



everywhere in this section we introduce: 

e = 4 (20) 
for Burgers, while e = 4 — d for the manifold allows perturbativc control. The rescaled velocity correlation satisfies 



td t A a p{\i) = 

Using A a/3 (u) = -d a dpR(\i) 
of the particle problem 



A a p(u) - ^d a d fj A 7<5 (u) - A jS (0) 



(21) 



-ii^g(u), it can be recast as an equation for the correlator of the random potential 



1 



td t R(u) = - - 2C + £u • V„ R(u) + - R'UuY - 2RU0)R'Uu 



(22) 



One recovers the familiar zero-tempcrature FRG equation for manifolds, derived here directly for the inviscid 
Burgers problem in N dimension, i.e. for d = 0, by identifying tdt = — \md m (for N = 1 it was obtained to 4 loops 
in [61, 62]). 

Let us focus on isotropic turbulence and denote i?(u) = h(u), the general case being very similar. The 
first property of the above FRG equations is that as long as the velocity flow is smooth, A 7 ^(u) — A 7 ^(0) ~ 



'■^h""(0)(u 2 S 1 s + 2u^us), hence from equation (18) the energy is conserved, 
d t £(t) = ^d t A aa (0) = 



(23) 



in agreement with standard knowledge for Burgers flows. However it is known (since Larkin, for review see e.g. 
[90] in the context of elastic manifolds) that h""(0) diverges in a finite time t c . Furthermore it is known since 
[57, 77, 91] that the solution of the one-loop equation At )Q| g(u) develops a cusp at the origin, more precisely 
h(u) = h(Q) + ^"(O) 2 !- + ^'"(O) 2 !- + a property which was found to hold also to next order (two loop) [89] and, 
from the physics of shocks, is believed to hold to any order. Hence in the present context of Burgers turbulence it 
implies non-conservation, and decay, of the kinetic energy, 

1 



W) = - dt A aai o) = -t.^L±V h »W - -ffi±%"(o)V- 



(24) 



If Note that we have assumed that A a f,.j,(c5) — Ai, a .;,((5) is continuous and vanishes at S = 0. 
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using that now A 7 a(u) — A 7 a(0) ~ — ^h'" {Q)u{8 1 s + u^us) where in the last equation we have substituted the 
scaling form (19) valid at large t. In the large-time regime, At i0l ^(u) flows to a fixed-point function A*^(u), which 
represents an asymptotic self-similar decaying solution. The energy decay can also be written as 

£(t) = ±A* aa (0)t^ , (25) 

an exact relation for the amplitude if one knows the fixed point to all orders. To one-loop accuracy this is equivalent 
to (24). Equation (25) generalizes the result (324), or (325), in [62] to any dimension N. 

As is well known from studies of the FRG equations [57, 77, 91], there are two types of fixed points. First a 
family of LR fixed points, such that at large u, A Q ^(w) ~ l/u n+1 , and R(u) ~ u n can be obtained by neglecting 
the non-linear terms (which are subdominant at large u) and considering only the linear part (rescaling) of the 
FRG equation (21). This easily recovers £ = Clr(^) = e/(3 + n), the result discussed in Section 3, and the so-called 
persistence of large eddies. Second, the SR fixed point, for which only one value of £ is possible, and which is 
obtained by shooting in the fixed-point equation (21), (22) from u = 0, asking for a fast decay of R(u) at infinity. 
This gives a non trivial ( = Csr = jrjy + ^ N > wncrc ^ N decreases exponentially at large N [77] . The LR behavior 
holds up to £(n c ) = C SR; hence for small e and large N it suggests n c ~ N + 1. 

On the other hand, we know that for d — and any N, the analog of the SR fixed point should be the one 
given by Kida, i.e. £ = 1, and (up to the Gurbatov LR-SR crossover) that n c = 1 separates LR from SR. Hence, we 
see that while the LR regime is well captured by the FRG, i.e. the loop expansion in powers of e, the SR exponent 
C sr and the SR FP of decaying Burgers is not well approximated. One reason seems to be that Kida physics is 
controlled by rare events and extremal statistics, and seems to be better captured by the replica-symmetry breaking 
(RSB) method, which even leads to some exact results for the Kida FP [74, 75] (these can be extended to any N). 
In fact even the n = 1 marginal case also involves some replica-symmetry breaking physics (as well as a connection 
to random matrix theory) [75] . 

Of course the above discussion concerned scales larger than £(t). For u < t^ 2 , i.e. in the inertial range, the 
FRG gives the correct physics of shocks and energy transfer, with a cusp in A. To which extent the agreement 
(shock size distributions, etc.) can be made quantitative remains to be worked out in detail. 



4-3. FRG equation for Navier-Stokes in momentum space 



For Navier-Stokes the one-loop beta function is non-local in real space and thus easier to display in Fourier space 
(for a real-space expression see Appendix A). While the general case is displayed in Appendix B, we give here an 
expression valid for the subspace of flows such that A Q/ g(k) = P?g (k)A(k). For later use we introduce the potential 
R(k) such that A(k) = k 2 R(k) In real space we can also write A a p(u) = Pj /3 (9 u )A(u) = -(<5 Q/ 3V 2 — d™dP)R(u). 
The function A(k) is then the Fourier transform of the trace A aa (u)/(iV — 1) = A(u), and the potential function 
R(k) is the Fourier transform of R(u). 

For N — 2 all incompressible tensors can be written in this form, and this is not a restriction; we can even 
use discrete Fourier sums. For N > 2, this requires A(k) = A(fc) i.e. isotropic turbulence, and the k-continuum 
limit, i.e. an infinite box; the sums below thus become momentum integrals, as explained in Appendix B. The FRG 
equations are 

9 t A(k) =^- T ^6 k , k _ q , q [A(q)A(k-q)-A(q)A(k)] 1 (26) 
q 

= k lw^$ { {k2 ~ q2) [ (k - q)2 - ^ + {N - 2)fc2(k - q)2 l ■ (27) 

Note that 6k,k- q , q := —P c jm(k)Pjbc(k — cD-Pmb^) i s not invariant under q-)k-q, hence the first term can also 
be written in a symmetrized form, given in Appendix B. We now use the rescaled variables 

A t , Q/3 (k)=^- 2+Ar iA tiQ/3 (k^ 2 ) (28) 

to obtain 

rC C, „ \ , 2 



t&A(k) = f2-C-iV|-|k-V k J A(k) + j^-j J2 ^,k-q,q |A(q)A(k - q) - A(q)A(k)J . (29) 
^ ' q 

Let us point out that this one- loop FRG equation, i.e. the unrescaled form (26), is very similar to the so-called 
Quasi- Normal approximation (QN). For N = 3, one can check that one recovers here the direct v = limit of 
equation (VII-2-9) of Ref. [6]. Let us however point out that the spirit here is a bit different. First, we are looking 
at v = directly. Second, we are searching for a fixed point valid for all k, with the appropriate choice of £. Third, 
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we use these equations as a first step in a systematic renormalizcd small-time (i.e. loop) expansion, which must be 
analyzed before carrying out the loop-expansion program. In addition we have kept N, the dimension of space, 
arbitrary. 

In the turbulence "closure" literature one often sees quoted the EDQNM, the Eddy damped quasi normal 
approximation, which is believed to improve on the QN. It amounts to replacing in the above FRG equation 
^k.k-q.q &k,k-q,q^fc.p.g- where Ok,p.q = ^/t^k,p,q are the "eddy damping rates", phenomcnological parameters of 
the theory, a standard choice being Hk,p,q — Mfc + + f-q? there are two choices for /i^ either [6] [i 2 = k 3 E(k) or 
fj, 2 . = cii(Jq dpp 2 E(p)) 1 / 2 . An interesting question, left for the future, is to understand how the next order in the 
systematic (renormalized) small-time expansion would compare with these phenomenological extensions. 

4-4- FRG equation for Navier-Stokes in real space: isotropic turbulence 

The flow equation (29) for an isotropic solution can be rewritten in real space. As for Burgers wc introduce the 
rescaled correlators via 

R(u) = t^-iRiut-^ 2 ) , A(u) = ^- 2 A(ut" c / 2 ) . (30) 

For isotropic turbulence, R(u) = R(u) and A(u) = A(u) (no numerical factor is introduced for NS). The flow 
equation is simpler in terms of the function i?(u), since a large part of equation (29) is almost local as a function 
of i?(u). We parameterize the solution as 

R(u) = r(y = u 2 /2) , A(u) = -V 2 a R(u) = - [Nr'(y) + 2yr"(yj] . (31) 

Then equation (29) turns into 

td t A(u) = (2 - C)A(u) + ■ 9 u A(u) + <5A L (u) + <5A NL (u). (32) 

The Laplacian of the first nonlinear term in equation (29) is local in real space and reads 

V 2 5A L (u) = 2 (3 -TV) N(2 + N) 2 r"(yf -8y 2 r {3) (y) (57V 2 + 37V - 40) r (3) (y) + (8JV - 20) y r (4) (y) 

-4yr"(y) (2 + N) {IN 2 + 3JV - 44) r (3) (j/)+4y ( (47V 2 + 67V - 26) r (4) (y) + (3JV - 5) y r (5) (y) ) 

- 2 (AT - 1) r'{y) (4 + N) {{2 + N) (jVr (3) (y) + 6 y r (4) (y)) + I2y 2 r (5) (7/)) + 8 y 3 r (6) (y) . (33) 

The second nonlinear term in the sum in equation (29) is strongly non-local in real space. Performing the angular 
average we obtain 

£A NL (u) = - e tk »R(k) A(k,p)p 2 R(p), (34) 

k P 

where 

fc 4 (4 + jV(4A-9))-4fc 2 (jV-l)p 2 + jVp 4 (fc 2 ~P 2 ) 3 „ (I 1 N 4k 2 p 2 \ 

A ^ P) = 2N{N-l) + 2(N-l)(k 2 +p 2 ) 2Fl {r 1 'Y>(k 2 +p 2 ) 2 ) ■ (35) 

This expression considerably simplifies for N = 3, with 

(I 3 4fc 2 p 2 \ fc 2 +p 2 / 2kp \ 

2Fl U' ' 2' (k 2 +p 2 ) 2 ) = -^kp~ atanh (PT? ) ' (36) 
and especially for N = 2 with 

see equation (67) below. 

4-5. Energy conservation and energy anomaly. 

Let us note some properties of the FRG equation for NS. First, as long as A a ^(u) is analytic at u = one has 
dtA. a p(0) = 0, which implies energy conservation. It can be seen by integrating equation (26) over k, and relabeling 
p = k q in the first integral and p = k in the second (also changing q > — q there), that 

* e A ( k ) =^tE bp+q r 2( : 6 + p y' q A(q)A( P ) = , 

k p,q 
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since the integrand vanishes by symmetrization of p, q. Since each integral contains terms of the form 
q Ql ...q Qi A(q), with i = 0, ...,4 (38) holds if A Q , l g(u) is smooth enough. 

Let us now recall where the energy anomaly comes from. The exact equation for n = 2 in the NS hierarchy 
implies the exact relation 



\dt{v« t ■ v ut ) = ,V^(v u , • v ut ) + jK,<«, 4 - < t )(v u , t - v ut ) 2 ) . 



(39) 



In the limit v — > it can be read either in the dissipative region u' — u (where the velocity is sufficiently smooth 
and the cubic term can be set to zero); or in the inertial range with u' — > u, where the first term is negligible. This 
expresses the energy decay rate as 



d t £ = -e = lim ^(v ut • V 2 v ut ) = Urn -V£ <«- t - < t )(v u , t - v ut ) 2 ) , 

v— >0 u'->u 4 

with £ = ±A QQ (w = 0). The relation (39) can be derived from equation (A. 2) in the Appendix which leads to 
<9 t A QQ (u) = -2V2^(u,u,0) = -2V» ut • v ot ) 

noting that V"/(v^ t v u / t • v U ' t ) = using translational invariance and incompressibility and that vanishes at 
coinciding points. 



(40) 



(41) 



Hence if there is an energy anomaly e > 0, the above implies C^^(u, u, 0) 



jjrVLy at small u. In the case of 



isotropic turbulence, using incompressibility, the third-order tensor can be parameterized by a single function of 
the distance h 3 (u), with h 3 (0) = as [79, 92] 



h 3 (u) 



l Q 0/3 7 



1 r 

2 



(A 



^.h 3 (u) ,.,.], . . . uh'Ju) — h 3 (u) 

1) \-h 3 {u) [Up5 a7 + U 7 5 Q(3 ) + U Q U,gU 7 2 



In the small-distance limit 

,(3) 



N, 



2c 



C^ 7 (u, 0, 0) ~ ^(0)u Q ^ 7 - Y ^(0)(u^ tt7 + u,5 af3 ) , ^(0) = N{N + 2){N _ 1) 
This is often expressed as 



-4e 



N(N + 2) 

and in particular leads to Kolmogorov's 4/5 law (for N = 3), 



42e 



N(N + 2) 



Note that for isotropic turbulence the two-point velocity-correlation can be written [79] as 



A a /j(u)= /(«) + 



uf'{u) 



N - 1 



Sal: 



uf'iu) ft 



-u a u fj 



(42) 
(43) 
(44) 

(45) 

(46) 



N - 1 

where the function f(u) defined in [79] is related to the potential function R(u) = h(u), as f(u) = —(A— l)h'(u)/u. 
Using equation (A. 4), one obtains the exact relation between the flow of f(u), and the third-moment function h 3 {u) 
[79], 



dtf(u) = - 



(N-l)h' 3 (u) + 



N 2 — 1 



h 3 (u) 



(47) 



generalized to any A, which recovers the above value for e using d t £ = ^d t f(0) = — ^N(N — 1)(A + 2)^,3(0) = — e. 
Now, in principle, from equation (A. 13) we have an expression for the third-order tensor (42), to lowest order in 
the (rcnormalized) small-time expansion, hence we can in principle relate the function h 3 (u) to A(u) (within our 
lowest-order FRG). The expression however is highly non-local and the dissipation rate ~ ^3(0) is not easy to 
calculate in general (while in Burgers it is, to one loop, simply proportional to A'(0) 2 ). 



4-.6. Enstrophy conservation and enstrophy anomaly in dimension 2. 

As is well known in dimension 2, A = 2, the velocity field is sufficiently regular so that the energy is conserved (no 
dissipative anomaly), and the energy flows to large scales (inverse cascade). There is no energy cascade towards 
small scale. We will show below that the solutions of the FRG equation satisfy these properties. 

In dimension A = 2 one also considers the vorticity field uj ut = e Q /39 Q v^ t . Taking the curl of the unforced 
Navier-Stokes equation (2), one gets 

9 t w u t + v ut • V u w ut = vV„w u t - /ACJ ut . (48) 
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We have temporarily added a friction term /i, which is often used to model dissipation at large scales. This implies 
that 

d t f f(u u t)= f f'(uut){vKuut-l^ u t) , (49) 

J u J u 

since the convection term integrates to a surface term, using compressibility, which is discarded. Hence there is 
conservation of any power, or function, of the local vorticity in the limit v — > 0, provided the right-hand-side has a 
vanishing limit. In particular 

(d t + 2m)(u&) = 2i/(w ut V»««t) , (50) 

hence under regularity conditions in the inviscid limit (existence of R^(0) is sufficient) one finds that the enstrophy 
D(t) = 5(^0*) i s conserved (setting friction to zero), 

d t D{t) = -\d t Vlm = . (51) 

Note that (50) is the limit v! — >• u of the more general relation for the decaying Navier-Stokes equation 

(dt + 2/i)<w u tW u / t ) = 2vVl (w ut w u / t } + ^VS((vS t - v2, t )(w», t - w ut ) 2 ) , (52) 

which allows to relate enstrophy non-conservation (the enstrophy anomaly) to the non-smoothness of the flow. The 
argument leading to (49) can be generalized to show that 

d t (e x ^) = —/i\d\(e Xul " t ) + v^lu^e^) . (53) 

In the absence of an enstrophy anomaly the limit v — > of the last term vanishes, and the solution is 

Z t (X) = Z t=0 (\e-^) ; (54) 

hence for /i = the full probability distribution of local vorticities is conserved. 

Until now these relations were exact. One now checks that the FRG equation does conserve enstrophy; indeed 
for a sufficiently smooth A(u) one has 

<9 t ^fc 2 A(k) = . (55) 

k 

This is shown as above from Sym p q [(p + q) 2 & p + q . p . q — _p 2 fr p , P + q .- q ] = 0, and is valid for N — 2 only. 



5. Short-distance analysis, cusp or no cusp? 

Here we study the behaviour of the velocity correlation A a ^(it) — A Q , ( g(0) ~ |u|'» 2 at small u, i.e. in the inertial 
range predicted by the FRG equation. To prepare for Navier Stokes, we first analyze Burgers, where we already 
know that the singularity is a linear cusp, i.e. £2 = 1, and we study it in Fourier space, since Navier Stokes is easier 
to express in Fourier space. Note that the analysis below can only exclude a range of values of £2, but to confirm 
that the selected values do occur in the solution, one must solve the fixed-point equation. 



5.1. Burgers 

For decaying Burgers, the FRG equation reads, in Fourier space 

^A(k) = (2 - C ~ n£- C -k ■ V k ) A(k) + £ fc2 2 [ a T f [ k J a ^ )]2 A(q)A(k -q)-£ ^-^A(q)A(k) . (56) 

It is true for any N and any symmetry (periodic sums, etc.). Note that while this equation seems to be non-local 
in real space, if one expresses it first using A Q ^(k) = P^(k)A(k), then performs the Fourier transform, it becomes 
local as a function of A Q( g(u). 

We want to understand why there is a linear cusp, and why there can be nothing but a linear cusp. For that 
we start with the isotropic case and look for a solution which at large k takes the form 

A(k)~AG(fc) , G(k) =k- N -^g(m 2 /k 2 ) (57) 

with g(0) = 1, and in Appendix F we note b = N + £ 2 - A further assumption is that g(p) admits an expansion of 
the form 1 - ^i% 2 + the mass m parameterizes the amplitude of the leading subdominant term, and its value 
is unimportant for the following + . A convenient heuristic form in that case is G(k) = (k 2 + m 2s j-( N +^)/ 2 . One 

+ Functions g(p) = 1 +p a can also be tried but they lead to additional conditions hence do not change the main point of the discussion 
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first notes that (57) implies that for (2 > 0, both integrals in (56) are UV convergent, thus A(u = 0) exists. The 
cusp corresponds to ($ = 1- In principle we can start by restricting our search to C2 < 2, i.e A(u) has no second 
derivative. If one tries to expand the first part of the non-linear term in equation (56) for q <C fc, one finds that 
it cancels half of the second term, and that corrections are of the order of J2 q 9 2 A(<z), which is a divergent series. 
The second half of the second term comes from the region q — k <C k. Hence the integral in (56) is dominated by 
large q. One can thus insert the ansatz (57) into equation (56), which shows that, up to the global factor of A 2 , 
the non-linear term behaves as 

f V\ ( k- q) f G(q)G(\k q|) - ^t G (q)G(k) = Bk 2 ^- N + Ck 2 - 2 ^ N + Dk^ N + ... (58) 
Jq 2 q (k q) q 

It is easy to see that B cancels between the two non-linear terms, see Appendix F.l, and that for the ansatz (57) to 
be a solution for Q 2 < 2, we need that the amplitude C = C(7V, £ 2 ) vanishes. The latter is computed in Appendix 
F.l, equation (F.6) as 



; [jv + C 2 (2 + c 2 )]r(-%)r(f + c 
r(i^)r(2±^±c,) 2 

It has a unique solution £2 = 1 in the interval £2 G [0, 2], and thus there only is a solution with a cusp. 



c = -i-jv+c 2 „.i/2-jv/2 l ly S2y^ -r V2;j * y - -2)i y T V*) = q / 5g % 

r( 1 ~ C2 )r( 2+jV+C2 f ' 



5.2. Navier- Stokes 

Inserting the ansatz (57) in the Navier-Stokes FRG equation (29), since ^k.k-q.q ~ k ■ q at large q, the second 
integral is UV convergent for C2 > 0. An analysis similar to the previous section yields 

]T 6 k)k _ q , q [G(g)G(|k - q|) - G(q)G(k)] = Bk 2 ~^- N + Ck 2 ~ 2 ^- N + Dk-^- N + ... (60) 
q 

The contributions to the amplitude B again cancel between the two integrals, hence B = 0. The amplitude 
C = C(N, (2) is computed in Appendix F.2: 

2fcpr - 2)c 2 - (N - i)](N + c 2 )r(-%)r(f + c 2 ) 4^F7vr(f ) 



c 



4(4 7 r) Ar / 2 



r(^)r(2±^)2 sin (i!|a)r(i±^)r(^) 



. (61) 



Solutions of this equation are plotted on figure 1. First, one can check that for any N, the function C(N, (2) 
vanishes for (2 = 1: C(N, 1) = 0. Next one finds that for N > A2 = 2.1155, the root C2 = 1 remains unique in [0, 2] 
and the linear cusp is the only possible solution. For N < N2, an additional pair of solutions appears on both sides 
of C2 = 0.87.... The largest one reaches 1 at N3 = 2.1145, while the smallest root is £2 w 0.8. For N < N3 the two 
additional roots are on both sides of 1 and, as TV — >• 2 + , one reaches while the other reaches 2. For TV = 2 (and 
N < 2) C(N, (2) is decreasing as a function of C2 for (2 € [0, 2] and the cusp is again the unique root jj. 

In addition to the cusp, there are other roots with £2 > 2. For N > 2.43... one finds that C vanishes exactly 
once in each interval (2 G [2p, 2p + 2], p = 0, 1,2, ... and diverges to ±00 for (2 = {2p) T . In the other intervals, 
p > 1, the root tends to ( 2 = 2p + 1 for large N. For N = 3 the other roots are at C2 = 3.32358, 4.98205, 7., .... For 
2 < N < N4 = 2.43... the root for (2 G [2, 4] may not exist in which case there is a double root in the interval [4, 6]. 
One root crosses ( 2 = 4 at N = 7V 4 . For N = 2 the roots are C2 = 5.02421, 7.0006, .... 

A peculiar result is that as N — > 2 + one root tends to 2~, but then seemingly disappears for N = 2, suggesting 
that this case has to be treated with more care. Indeed, we will see in section 6.2.4 that N = 2 and £2 = 2 is indeed 
a solution. The calculations are rather non-trivial, since the integrals are not defined without proper regularization. 

To conclude on an optimistic note, although the cusp seems to be the only solution for TV = 3, we did find 
some non-trivial values for Q 2 for N slightly larger than 2. One possible scenario may be that these become valid 
in a larger domain in N when higher loop corrections (higher powers of time) are included. 



6. Two-dimensional decaying turbulence (N = 2) 

6.1. Basic properties 

In dimension N = 2, since Kraichnan and Batchelor it is believed that [41, 6]: 

* Note however that there are other roots, i.e. £2 = 2p+ 1, p = 1, 2, ... which arc potentially possible behavior. For N = 1 such solutions 
are A(it) — A(0) = 5Zm=i £t2m^ 2m + i2p+lM 2p+1 + and are formally possible solutions at small u, but do not seem to correspond 
to globally physical fixed points. 

J For N < 0, there are again additional solutions for £2 6 [0,2], but for NS these roots are not physically interesting. 
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Figure 1. Blue solid curves: Locations in the (C2;^V) plane, where C given in equation (61) vanishes. Orange 
dashed lines: N = 2 and N = 3. For N = 3 the solutions in [0,4] arc Ca = 1 and ( 2 = 3.32358. For N = 2 they 
are £2 = 0, 1, 2. Note that the maximum of the lowest curve is at N = N2 = 2.1155, £2 = 0.87... hence in the (very 
small) interval 2 < TV < 7V2 there are two roots on both sides of this maximum (in addition to the root at £2 = !)■ 
At N = N3 = 2.1145 the upper root reaches £2 = !■ 

(i) the energy flow is to small k while the enstrophy flow is to large k; 

(ii) since there is no direct energy cascade there is no energy anomaly, i.e. lim„_ > o+ fA"(0) = 0, and energy is 
conserved i.e. <9 t A(0) = ±d t J dkE{k,t) = 0, once one neglects dissipation on large scales due to e.g. friction; 

(iii) there is an enstrophy cascade, i.e. there is an enstrophy anomaly lim„_ s . + vA""(0) ^ and enstrophy is 
not conserved -«9 t A"(0) = ^d t J dkk 2 E(k,t) < 0. 

Let us see how these feature arise from the FRG equation. To facilitate the calculations, we introduce the 
stream function ip and the vorticity u such that v^ t = e a 0dpip u t and u) u t = e Q ^9 Q v^ t = — V 2- u t. In terms of the 
stream function, R and A aj g are 

R(u - u') = (VWu't) , (62) 
A a/3 (u) = -(<WV 2 - V Q V /3 )i?(u) . (63) 
For isotropic turbulence A(u) = A aa (u)/(N - 1) = A QQ (u) and A(fc) = k 2 R(k) = E(k)/k. 



6.2. Isotropic turbulence 

6.2.1. FRG equations. In section 4.4 the FRG equation has been given in real space. We remind that one 
parameterizes the isotropic solution as 

R(u) = r(y = u 2 /2) , A(u) = -V 2 u R(u) = -2 [r'(y) + yr"(y)] . (64) 

For isotropic turbulence and N = 2, equations (33)-(35) simplify to 

td t A(u) - (2 - C)A(u) + CydyA(u) + SA L (u) + 5A Nh (u), (65) 



6A L {u) = 8 
5A NL (u) = 



r"{xY dx - 2r\y) r"(y) + 3yr"(y) z - Ay r'(y) r (3) (y) + y 2 r"(y)r (J, (y) - y 2 r'(y)r w (y) 



Mi 



Mi 



-ik-u 



R(k) 



{k 2 - q 2 ) 2 q 2 R{ q ) 



J\<l\<k 

k qdq 



, (66) 
(67) 



where Jiq^j. '■— J %7- The integration constant has been fixed assuming that A(u), and thus the above correction, 
vanishes at |u| = oo. 

The above can be turned into an equation for r'(y): 

2r'(y)r"(y) - yr"{y) 2 + yr'(y)r'"(y) - 3 f dx r"(x) 2 ^ 



td t r'(y) = (2 - C)r'(y) + (yd y r'(y) + A 



A d f 1 f°° f 2 ^ 

+ -— / dA (A — 1) / dz d9 
I" dy Jo Jo Jo 

f(x) = [xr'(x)}' = r'(x) + xr" (x) . 



f{y + \z + 2^MTz cos 6) - f(y) f(z), 



(68) 
(69) 
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In Fourier space, the FRG equation reads 

td t A(k) = (2 - 2C)A(fc) - |fcA'(fc) + SA(k) . (70) 

Using the distance-geometry representation of Appendix E, the sum of the two non-linear terms SA(k) = 
S Al (fc) + <5 Anl (fc) can be rewritten in Fourier space as 

1,4 poo pi , j, r t, - 

,5A(fc):= — J ds J ^ dt [(s tf - 4] -_ V(« 2 - 1)(1 - t 2 ) A( - (s - t)) A(fc) - A[-(s + *)) I . (71) 
We now study the scaling form and the asymptotic behavior of the fixed-point solution A* (fc) for large and small fc. 

6.2.2. Searching for an isotropic fixed point We look for a fixed point of the form 

A{u) = t c - 2 A{ut-^ 2 ) , A(fc) = f 2C - 2 A(fci c / 2 ) . (72) 

The asymptotic behaviors at small and large distances are 

A(0) - A(u) ~ u C2 for u < 1 <=> A(fc) w Afc~ (2+C2) for fc > 1 (73) 

A{u) ~ for u > 1 A(fc) w ^fc™" 1 for fc < 1 . (74) 

valid only up to logs (and for A{u) at large u only an upper bound since for integer n a faster decay is possible 
from analyticity in Fourier space). 

Consider now the mean kinetic energy E(t) = 27r(v 2 ) which is given in the non-standard units introduced in 
equation (8). If we assume scaling, then 

poo poo 

E(t) = A(k)kdk = t 2 - C A( P )pdp . (75) 
Jo Jo 

If C2 is large enough, the energy should be conserved (the cusp seems necessary for the violation so let us consider 
C2 > 1)- Then the value naturally compatible with energy conservation is £ = 2. In the context of disordered 
system this is called the Larkin exponent, i.e. the dimensional reduction exponent £ = e/2 with e = 4. 

Furthermore Batchelor [41] proposed that E(k,t) — v 3 tf(kvt) which implies conservation of energy, if J dzf(z) 
converges. This is again ( = 2. Let us recall that E(k,t) = fcA t (fc). It also implies a decay of the total enstrophy, 
i.e. of / dfcfc 3 A t (fc), proportional to t~ 2 if J dzz 2 f(z) converges. Assuming that E(k) is independent of v at large 
fc implies f(x) ~ x~ 3 and E(k,t) ~ t~ 2 k~ 3 at large fc. This leads to 

C2 = 2 . (76) 

This is indeed the only solution we found to be consistent with our analysis of the flow equations ft as we discuss 
below. At small fc it behaves as A t (fc) ~ t A k 2 . Note that it is not a long-range fixed point with n = —1, which 
would also give £ = 2 according to the general discussion of section 3.2. The reason is that the amplitude of the fc 2 
depends on time (in addition the energy would be infinite). Rather it corresponds to n = 3, but a short-range fixed 
point. 

We now consider the general properties of the small- and large-fc expansions of the fixed-point solution A*(fc). 
Since the fixed-point equation is neither local in real space nor in Fourier space, the expansions of interest are 
expected to contain unavoidably global properties of the fixed point A*(fc). 

6.2.3. Small-k expansion. The expansion of the nonlinear term (71) in the distance geometry representation is 
given in Appendix G. To lowest orders it reads 

fc 2 f°° fc 4 f°° 

<5A(fc) = — / dqqA(q) 2 -— dqqA'( q ) 2 + 0(k 6 ) 7 (77) 



47T 7o 167T 







where we used the small-fc behavior (74) that implies, for n > 2, A(fc = 0) = A'(fc = 0) = 0. Assuming that 
A*(fc) = A-2k 2 + A^k 4 + ... we find from the FRG equations in Fourier space: 

1 r°° 

(3C - 2)A 2 =4A 2 = — dfc fcA(fc) 2 , (78) 
4tt Jo 

1 f°° 

(4C - 2)A 4 = 6A 4 = — — / dfcfcA'(fc) 2 . (79) 
167T J q 

ff from the FRG equation in Fourier at large k one sees that S A* (k)/ A* (k) should go to a constant at large k equal to — (2 — ( + 5 CC2)- 
This can be checked numerically and we found it holds only for £2 = 1,2 and that £ 2 = 1 can be excluded as it would be compatible 
only with ( > 4 a value much too large 
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ln(jfc) 



Figure 2. Double logarithmic plot for the rescaled energy E(k) := kA(k) as a function of k. The solid black line is 
our solution of equation (65); the thick dashed lines show the asymptotic slopes ±3. The thin grey-blue dashed line 
is Bguess(fc) := fcAguess(fc), see equation (85). 



6.2.4. Large-k expansion. As shown in Section 6.2.2, the large-fc, i.e. small-u asymptotics of the fixed-point solution 
is given by equation (73). Assuming that A(fc) = A/k A + <3(l/fc 6 ) we find in Appendix H, that the nonlinearity 
gives 

+ 0$*). (80) 

Inserting equation (80) into the flow equation with £ = 2 and asking that it be at a fixed point yields 

A = 16tt o A(fc)~^ (81) 

fc 4 

at large fc, which implies A(fc) ~ Wcz, noting that this power of fc is preserved by rescaling by t^l 2 in (72). This 
means that in real space, using the definition (31) to pass to the second line: 

A(u) - A(0) = 2u 2 lnu + ... (82) 

R(u) = -A(0)^-iu 4 lnu+... (83) 
4 8 

6.2.5. Numerical solution To find numerically a fixed-point solution A* (it) of equation (65) is highly non-trivial. 
Since equation (65) is an integral-equation, none of the standard techniques, such as Taylor-expansion, or solution 
as an eigenvalue problem are available. From a decent physical fixed point, we expect that it is attractive w.r.t. 
all (sensible and small) perturbations. Thus, if we propose a guess A guoss (u), which satisfies the above mentioned 
asymptotic behaviors and constraints, and is close to the true solution, it should converge against the fixed point. 
This is what we succeeded in doing, starting from A cor (fc) = 0: 

A(fc) = A gucss (fc) + A cor (fc), (84) 

fc 2 54.7237fc 6 2.65177fc 4 ~ 

( (fc 2 + l) 3 + (fc 2 + l) 6 + (fc 2 + l) 5 

The main problem then was that numerically the flow-equation (65) is rather unstable. The technique we finally 
succeeded in getting to work was: Starting with A corr (fc) = 0, we recursively inject A corr (fc) into the flow equation 
(65), and use the latter to evolve A(fc) during a small time-step, giving us an improved approximation for A(fc), 
calculated for approximately 100 fc-points. The latter is then projected onto an optimal spline with only 20 
supporting points, or more precisely a non-linear transformation thereof. This procedure is numerically much more 
stable than using a best polynomial fit, a Fourier representation, or any of the other known sets of orthogonal 
functions we tried. The projection effectively smoothes the function, while still capturing the necessary details. 
The complete technical details can be found in Appendix I, most importantly a check of the convergence of the 
function, see figure 12, as well as its tabulated values. Here we illustrate the result in the form of the classical 
double-logarithmic plot for the energy as a function of fc, see figure 2. The small-fc asymptotics is E(k) ~ fc 3 , 



A guess (fc) = 16tt ( — + — + — . (85) 
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and the large- fc asymptotics is E(k) ~ k 3 . We remark that the solution remains below the asymptotic small- k 
behavior, but then converges from above towards the asymptotic large-fc behavior. 

6.2. 6. Physical interpretation of the solution We have found a fixed point A(fc) for 2D decaying turbulence. Having 
in mind equation (72), the velocity-velocity correlation is A(fc) = t 2 A(tk) where A is time independent. Thus 

td t A(k) = 2A(k) + fcA'(fc). (86) 

This implies similar relations for the time evolution of energy E(k) = kA(k) and enstrophy f2(fc) = k 3 A(k) 

td t E(k) = 2fcA(fc) + k 2 A'(k) , td t n(k) = 2fc 3 A(fc) + fc 4 A'(fc) . (87) 

We define scaled momentum-dependent energy and enstrophy decay rates, written in terms of the scaled momentum 
k = kt^' 2 = kt and correlator as: 

E(k) := d t E(k) = 2fcA(fc) + k 2 A'{k) = -d{j E {k) (88) 

tl(k) := t 2 d t n(k) = 2fc 3 A(fc) + fc 4 A'(fc) = -d{j n {k) (89) 
where we have defined the scaled energy and enstrophy fluxes: 

j E (k) = -k 2 A{k) , j n (fc) = -fc 4 A(fc) + 2 [ dpp 3 A(p) (90) 

Jo 

The scaled energy and enstrophy decay rates, and the scaled fluxes, all as functions of k are plotted on figure 3. 
(For convenience we revert to the notation of the rest of the paper for the argument of A, i.e. the a;-axis is 
called k but it is more properly k). We see that the energy flux jE{k) is negative, thus to small momentum 
scales, and moreover energy is conserved, since lim^oo is(fc) = 0. Hence the total energy decay rate vanishes, 

e = — dt£ = J dkE(k) = 0. The enstrophy flux is mostly positive, thus to large momentum scales, and 
enstrophy is not conserved. However, note that for A(k) ~ 1/fc 4 , the integral in (90) grows as lnfc, thus the 
enstrophy conservation is only marginally violated. For A(fc) ~ l/[fc 4 (lnfc) c ], and c > 1, the enstrophy would be 
conserved, since the flux at large k would vanish. Thus a small modification of the asymptotic behavior, which 
might not be given correctly by our leading-order fixed point, would be enough to ensure enstrophy conservation. 

To summarize, it seems that the FRG fixed point is compatible with the Batchclor-Kraichnan scenario [41, 93] 
with an enstrophy anomaly — dth(^ t ) = lim„^o ^((Vw) 2 ) = e u . Since e w has dimension (time) -3 , the energy 

spectrum within the Batchelor-Kraichnan 2D enstrophy cascade is E(k,t) ~ e^/ 3 ^ -3 , with, in decaying turbulence 
e w ~ 1/t 3 . However, it is also known that this is not the end of the story, and that more detailed arguments and 
assumptions lead to additional logarithmic corrections, e.g. E(k,t) ~ ef/ 3 /(fc 3 (ln k) c ) with c = 1/3 [93, 6]. It was 
also argued that in real space the vorticity correlations are not ((ui u — wo) 2 ) ~ lnw at small scale, as would be 
the case if E(k) ~ 1/fc 3 . Instead they have fractional powers of lnu, claimed to extend to all moments of the 
vorticity field, as ((oj u - cj ) n ) - (lnu)™/ 3 [94], as a consequence of the infinite number of conservation laws; each 
conservation law is violated and leads to a flux of the corresponding (almost conserved) quantity. (Equivalently 
one can write that (v ■ Vixj n uj n ) is a constant [94]). Most of these issues were discussed for forced turbulence, but 
remain relevant for the decaying case. Note that in forced 2D turbulence there is an additional regime with an 
inverse energy cascade E(k) ~ fc~ 5 / 3 (see e.g. [6]). The energy flows to large scales, until the largest scale is reached, 
where coherent structures form. (These lowest k modes may be called a condensate). At small scales however, the 
behavior should not be qualitatively different from decaying turbulence. For numerical and experimental tests of 
the enstrophy cascade see [95, 96, 97]. As discussed in [98, 99], friction is a marginal perturbation, hence should 
change the logarithms of the distance in the vorticity correlations into power laws. Finally, for a more mathematical 
discussion of the enstrophy anomaly sec [100, 101]. In particular, the anomaly was proven to vanish in the forced 
2D Euler equation with friction [102]. 

A challenging question is whether some of this physics can be captured in higher-loop extensions of the present 
approach. 



6.3. Periodic 2D-turbulence 

Let us consider the NS equation in a square box of size L = 2n with periodic boundary conditions. We study the 
FRG equation (26) for TV = 2, using integer Fourier modes (k x ,k y ) £ J?, which are summed over. It is easy to 
analyze numerically the FRG equation projected onto a grid [-Q, Q] 2 , setting A(fc) = outside. Equation (26) can 
be written schematically as d T R(k) — (R * R){k) with r = t 2 /2 and R(k) = A(k)/k 2 . Rcscaling is not crucial here, 
since periodic turbulence corresponds to £ = 0. One finds, for any Q, that the flow asymptotically behaves as 

R(k) = fi4=(±i,o) : (o,±i) + he~ Xr (91) 
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and that the energy becomes entirely concentrated in the lowest modes k 2 = 1. The Fourier coefficient r\ 
of these modes tends to a constant at large times, while all other modes decay. There is a transient regime 
where the other modes first increase before decaying, following (91). The are obtained by diagonalizing 
<9 T ffc = 2(Ri * r)(k) = — Affc, where Ri = 7 , i^fc=(±i,o),(o,±i) is a re& l non-symmetric matrix. We truncated 
numerically on a grid k € [— Q, Q] 2 and found, apart from one trivial eigenvalue A = 0, corresponding to R = R\, 
that all eigenvalues are negative. Only the leading one corresponds to a vector with all positive entries for k 2 > 1, 
which is requested from (91). A plot of A versus 1/Q 2 is approximately linear and the numerical solution suggests 
A = — 0.6ri. The leading eigenvector is very well fitted by fk/f\ — 1.57k~ a with a = 5.8 ± 1, which we checked up 
to Q = 16. Since a > 4, this is consistent with the absence of an energy anomaly. It is interesting that this value 
of a seems to be indeed near the Batchelor value of a = 6, which corresponds to A(/c) ~ k~ A and E(k) ~ fc~ 3 [41], 
consistent with our analysis in the last section. 

In conclusion, we want to note that R(k) written in real space tends to a fixed point which corresponds to the 
average over the set of exact, time-independent solutions of the Euler equation, 

iput = w\ x cos(u x ) + W2x sin(u 2 ,) + wi y cos(wy) + W2 V sin(u y ), (92) 
v ut = ( - wi y sin(uy) + w 2y cos(u y ), Wi x sin(u x ) - w 2x cos(u x )) . (93) 

This is easily checked by inserting into Euler's equation. The four parameters uii X . y are independent Gaussian 
random variables with zero mean and variance ri/2. The FRG suggests that the way it tends to this fixed point is 
non-trivial (with a non-analytic correlator) . It would be of great interest to study that question in detail. 

7. Analysis of the FRG equation in three dimensions (N = 3), and in large dimensions (N — > 00) 

Obtaining a numerical solution for the fixed point of the FRG equation for N > 2 is difficult. In three dimension 
N = 3 we have studied the FRG equation (26) for a periodic flow, in Fourier space, very much as we did in the 
previous Section for N = 2. We have found that on a grid in Fourier space [— Q, Q] 3 equation (29) does flow to 
a fixed point with £ = 0. This fixed point, to our numerical accuracy, was compatible with a cusp behavior at 
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large k, i.e. (2 = 1- Since the result is not too surprising, and consistent with the analysis of Section 5, we will not 
reproduce the details here. Instead we now turn to the large dimension limit (large TV limit) analysis of the FRG 
equation. 

Large dimensions are often a means of controlling an expansion. The aim would be to sum contributions from 
all loops at TV — > 00, thus avoiding any artifact from a closure scheme, as has been done for random manifolds, 
i.e., Burgers, [54, 103, 104, 105]. While this remains a project for future research, we have analyzed the one- loop 
equations for TV — > 00. In terms of the function r(y) introduced in equation (31), the RG equations at large TV are 
derived in Appendix J and reads 

TV (2 + TV) td t r"(y) + 4y [(2 + TV) td t r'"(y) + ytd t r""(y)\ = 2r"(y) + ( yr"'(y) + 2r"{yf + 2r"'(y) [r'(y) ~ r'(0)] 
~4r"(y) + y(8 + 6( + (- i y"(y) + 4( y 2 r""(y)+ 2r'\yf + 10r"'(y) [r'(y) - r'(0)] 



TV 



+ 12yr""(y)[r'(y) - r'(0)l + 28yr" (y)r'" (y) 



dtr'"(y + t)r"(t) + 0(1/TV 2 ), 



(94) 



where C — Co + C-i/N + ... We now look for a fixed point. To leading order we have 

= 2r"(y) + Coyr"'(y) + 2[r"(y)] 2 + 2r"'(y)[r'(y) - r'(0)], (95) 

which coincides, up to a numerical prcfactor, with the large-TV limit of the 1-loop Burgers equation, see e.g. 
equation (7.7) of Ref. [104]. This confirms that at least to one loop the infinite- TV limit of the decaying Navier 
Stokes equation reproduces that of the Burgers equation. Equation (95) has an exponentially decaying solution 
only for Co = 0; an analytic solution for the inverse function can be written as 

z := -r'(y), (96) 
y = z - z Q - z \n(z/z Q ). (97) 



The asymptotic behavior is z 



z e~ 



-v/zo f or l ar ge y. However, we cannot neglect the terms of order 1/TV for 



y/N 3> 1, and the above solution is valid only for y <C TV, the primary region. On the other hand, for y 3> z$, one 
can neglect the nonlinear terms in equation (95) due to the exponential decay of r{y). Both solutions are expected 
to match in the inner region zo ^ 2/ ^ TV, which becomes quite wide for TV — > 00 [77, 89]. Presumably in the inner 
region both solutions have a simple exponential behavior, to order 1/TV. The linearized equation to order 1/TV reads 

4r"(i/) + 8yr"'(y) + C- lV r'"{y) - 10r"'(y)r'(0) 



= 2r"(y) - 2r"'(y)r'(0) 



N 

12yr""(y)r'(0) 



dtr"'(y + t)r"(t) 



(98) 



We assume that this equation has a solution which behaves as e~ y ^ z ° for y ^ z . We are free to fix z = 1 for the 
sake of simplicity. Substituting r{y) = Ce~ v into the linearized equation, and collecting all terms proportional to 
e~ y and ye~ v , we obtain 

1 r 



ye 



= 2 




2/(0) 



N 



10r'(0) 



die-V'(t) 



1 r 

TV 



The first line gives r'(0) 
4 

TV " 



O 



C-i - 12r'(0) . 

1 + 0(1/ N) and the second C-i 
1 

TV2 



4. so that 



(99) 



(100) 



Close to Zo the solution of equation (97) has the form z 

c 2 = i. 



zq = \/2zQ~y ~ u. This implies that the cusp persists, i.e. 



8. Decaying surface quasi-geostrophic turbulence 

An interesting and still much studied generalization of 2D NS is the Surface Quasi Geostrophic (SQG) equation. It 
is defined in dimension TV = 2, and depends on a continuously varying parameter a. In real space it reads 

d t T ut + v ut • VT ut = vV 2 T ut , v ut = £ x V^ut , (-V 2 ) Q / 2 ^ut - T ut . (101) 

It describes the convection of a quantity T ut by the velocity field v ut , which in turn is related to the velocity. 
For a = 2 one recognizes that the quantity T ut is precisely the vorticity w u t, and one recovers the usual 2D NS 
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equation. For a = 1 the field T uf represents the temperature in the "true" SQG turbulence, which is used to model 
the 2D atmospheric flow on the surface of the Earth. Finally, for a = — 2 the model was obtained by Charney and 
Oboukhov for waves in rotating fluids, and by Hasegawa and Mima for drift waves in a magnetized plasma in the 
limit of a vanishing Rossby radius. It is thus called the Charney-Hascgawa-Mima equation [106]. The naive scaling 
dimension of the field T ut is 5 U T ~ u H with H = (2 — 2a)/3. Recently, it has been conjectured that isolines of T ut 
in the inverse-cascade regime of the forced Charney-Hascgawa-Mima equation were SLE lines with k = §(1 + 2a), 
with some numerical evidence [42, 107]. 

The SQG equation for arbitrary a shares some properties with the 2D NS equation, in the sense that in the 
inviscid limit both the enstrophy D = \{T 2 ) (and all powers of T) and the energy E = ^(Tip) are conserved for 
flows smooth enough. To show the latter one goes to Fourier space, where the relation between ip ut and v ut is 
V> q * = <T Q ?V This yields 



k,q 



q a k a 



q-kt^-qt 



o , 



(102) 



due to the symmetry q <R- — k. 

Here we consider the decaying inviscid SQG equation. We display the FRG equation to one loop, leaving its 
analysis for the future. One defines the 2-point correlation in Fourier space, 



A T (k) = (T_ kt T kt ) . 
The FRG equation for A-r(k) is derived in Appendix D, and reads 



(103) 



d t A T (k) =i^A T (q)A T (q + k)[qx k] 2 {< 



-2a 



2q~ 



k|~ 2a } 



-2^A T (q)A T (k)[qxk] 2 {q 
q 



<1 



— a i—a 



'k 



<r Q |q- 



Note that the correlator associated to the velocity is A(k) 
rescaled correlators via 



A T (k) =;- 2+ ^ 6 ~ 2a )A T (k^/ 2 ), 

A(k) =;- 2+2C A(k^/ 2 ). 



-° + £T a |q + kr a } . (104) 
A; 2-2a A<r(k). Again, it is convenient to introduce the 

(105) 
(106) 



For the rescaled correlator, the FRG equation can be written as 
i<9 t A T (k) = [2 - (3 - o)C]A T (k) - ^kd k A T (k) 



p~ a )A T (q) (q- a - p' a )A T (p) - 2(q- a - fc- a )A T (k) . (107) 



q,p=k-q 

One can use q x p = q x k. The equation for the rescaled velocity correlator reads 



td t A(k) = (2 - 2C)A(k) - |k9kA(k) 



E 



[qxk] 2 ~ 



A(q) 



(P a 



^ 2 A(p)-2 (pa - ga)(r ' ga) A(k) 



k a p a 



(108) 



q,p=k-q 

In this form it is easy to check that for a = 2 one recovers the FRG equation for the NS equation in N = 2. 

The FRG equation (108) written in real space for general value of a is a nonlocal intcgro-diffcrcntial equation. 
It is interesting that there are some values of a for which the FRG equation becomes quasi-local, i.e. involving only 
derivatives of finite order at the point u and at the origin u = 0. For instance that happens in the case of the 
Charney- Hasegawa-Mima turbulence corresponding to a = — 2. 

In Appendix F.3 we have studied, as we did for Burgers and NS, the possible values for the exponent C2, defined 
from A(k) ~ fc~( 2+l > 2 ) at large k, and isotropic turbulence. More work is necessary to study the fixed points of the 
FRG equation as a function of the parameter a. 



9. Conclusions and Perspectives 

In this article, we have applied functional-rcnormalization-group methods to decaying turbulence. In contrast to 
standard perturbativc RG, the functional RG approach takes into account a coupling function i.e. an infinity of 
couplings rather than one or few. It naturally leads to a non-analytic 2-point function. While the method is in 
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principle exact, as any RG treatment, in practice the flow has to be projected onto a lower-dimensional subspace, 
here the equal-time 2-point velocity correlation function. With this projection in mind, the FRG equations are 
organized in an expansion in powers of the 2-point velocity correlation function itself, equivalent to a loop expansion. 
As we have discussed, they correspond to a small-time renormalized perturbation theory. Here, we studied the 1- 
loop equations. For Burgers, they reproduce the FRG equations derived in the context of random manifolds, and 
correctly describe the singular structure of the flow, made out of shocks. While this had been worked out in details 
before for N = 1, here we extended it to any dimension N. 

Let us stress that the method works at least qualitatively for Burgers; that it correctly accounts for shocks, and 
that the distribution of velocities is not close to a Gaussian. The reason is that the extension to a manifold provides 
a model which can be controlled perturbatively (in d = 4 — e), while at the same time exhibiting shock singularities, 
non-conservation of energy (called failure of dimensional reduction in the context of disordered systems) and energy 
cascades. This is because shock sizes and the magnitude of the energy decay rate are O(e) in that expansion. That 
in itself is remarkable in the turbulence context, and motivated us to consider Navier-Stokes with this method. 

For Navier Stokes, the fixed point depends on the dimension. For N — > oo, the FRG equations converge 
(at leading 1-loop order) to those of the decaying Burgers equation. Thus the 2-point velocity correlation function 
should grow linearly with distance, i.e. have a cusp. This cusp is also the only possible solution for the 3-dimensional 
FRG equation, at 1-loop order, in contradiction to experimental evidence. It is possible, that at second (2-loop) 
or higher order, new non-trivial fixed points emerge. If this is not the case, one would have to understand why 
the method seemingly does not admit the correct singularities. Since for large N we find that the FRG equation 
reduces to the one of Burgers, hence has shock singularities, one possible way to understand that may be via a 
large-dimension expansion combined with a loop expansion. 

Finally, in two dimensions the FRG equations allow for a fixed point which is consistent with Batchclor's 
scaling. While the equations are similar to the quasi-normal Markovian approximation, we give here an explicit 
solution. Again it seems a good starting point to include higher-loop corrections, one challenge being to confirm, 
or infirm, the conjectured logarithmic corrections. We have also written the flow equations corresponding to SQG 
turbulence, which await a more detailed analysis. 

We hope that this work helps to bring a new perspective in a long-debated subject. 
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Appendix A. 1-loop FRG equation in real space 



The two equations from (13) needed to one loop in the inviscid limit are 



a t A QlQ2 (u) =Sym P ai ^(d u )C$ a2 (0, 0, u 



2Sym P Q T i/3 (a u )^C^ Q2 (0,0,u) 



Sym[dr4 3 > f0,0,u) 



-Sym 



Pa 1 ;^(5u)Cj 3 7 ) Q2 (u,U,0)l (A.l) 



(Euler) 
(Burgers) 



(A.2) 
(A.3) 
(A.4) 



with u = U12 and where P a .p 1 {d) = Pjp(d)d 7 + Pj 1 (d)d / 3 (Euler), P a -p y {d) = <5^ 7 <9 a (Burgers). In the first three 
lines Sym[...] means symmetrization over aa, a 2 and we have used that A QlQ2 (u) is even in u (no average helicity). 
In the last line, and everywhere below Sym[...] means symmetrization over u" 1 ,^ 2 ,^ 3 , i.e simultaneous exchange 
of the points in space and the indices. 

To lowest order we replace (denoting Uy := — Uj) 

C, Sa 2 c t3 ( U l' U l' U 2^3) = A /37 (0)A (u 23 ) + A Pa2 (u 12 ) A 7Q3 (u 13 ) + Ap a3 (u 13 )A 7Q2 (u 12 ) . (A. 5) 

Hence integrating equation (A.4) one gets 



C ^ a2a3 ( u u^2,u 3 ) = --tSym^P Ql/37 ((9 Ul )^A ( 3 Q2 (ui2)A 7Q 3(ui3)J + P ai ^(d Ul )^A Pa3 (u 13 )A ia2 {u 12 ) 



-3t Sym 



p J xl (dm ) (A/jc.3 (u 13 )<9^ A 7Q2 (u 12 ) + A /jQ2 (u 12 )<9^ A 7Ct3 (u 13 



-3tSym ^ i 1 (A ( 3 Ct2 (u 12 )A^ 3 (u 13 )) 



(Euler) (A.6) 
(Burgers) (A. 7) 
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where for Euler we used the transversality of A Q( g. Expanding, one finds for Burgers 

-f^CliUaa ( U l . U 2- U 3) = (A^ 2 (u 12 ) Ap a3 (u 13 )) + (A^ (u 21 ) A^ Q3 (u 23 )) + fl£ (u 32 ) A^ Ql (u 31 )) . 

(A.8) 

This expression is symmetric and does not need to be symmetrized. Taking the limit of u 2 — > Ui one finds 
-^'^^(u^u^ua) = d«l (A^ a2 (u 31 )A^ 1 (usi)) +A? a2 (0)d^Ap a3 (n 13 ) + A 0ai (0)fl£ A^ 3 (u 23 ) . (A.9) 
We have used that 

lim (aS 1 1 A (3Q2 (u 12 ))A^ Q3 (u 13 ) + (9S 2 2 (A^ 1 (u 21 ))A^ 3 (u 23 ) = , (A.10) 

U2— >Ui 

which comes from ^'^(u) being odd. One then gets 

a t A QlQ2 (u) = -td^d^A^(uf +2tA p,(O)A M , ai (u) , (A.ll) 

where we have used that A / g Q , 2;( g/ (u) is odd and that A / g a2;i g' ai (u) = — ^sfl' 0l03 (u) is symmetric in ai,a 2 . This 
gives the equation in the text. 

For Euler one finds by expanding 

+ Pl 27 (du 2 ) (A^ (u 23 )d£ a A 7Q1 (u 21 ) + A Pai (u 21 )d'i 2 A 7Q3 (u 23 )) 
+ ^a 37 (5u 3 )(A (9 a 1 (u 3 i)^ 3 A 7Q2 (u 32 ) + A^ 2 (u 32 )c^ 3 A 7Ql (u 31 )) . (A.12) 
This expression is symmetric and does not need to be symmetrized. We now take the limit u 2 > Ui : 
-^^(ux.ux.ua) = P Q T 2y (a U3 )(A^ /3 (u 31 )<A r7 (u 31 ) + A^(n 31 )d^ (n 31 )) 

+ u jim i ^ T y(3u 1 )(A^ Q2 (u 13 )a^Ay 7 (u 12 )+A /3 , 7 (u 12 )^A rQ2 (u 13 )) 

+ P 7 T 7 ,(a U2 )(A^ Q2 (u 23 )^A y(3 (u 21 ) + A w (u 21 )^ 2 A ya2 (u 13 ))J . (A.13) 

The term lim U2 _j. Ul (...) involves a non-trivial coinciding-point limit. One may naively equate it with 

Pj 7 ,(9 Ul )A^ 7 (0)<A 7 , Q2 (u 13 ) + P 7 r y (d u JA^(0)<A 7 , ct2 (u 13 ) , (A.14) 

but this is actually incorrect. It would lead to a term +2tAp> 7 (0)dP 9^A aiQ2 (u) in the beta function. The correct 
beta function must retain the non-trivial limit, for which we obtain, with u 3 i = u: 

a t A QlQ2 (u) = 2[Plp(d u )dZC$ a2 (o,o,u)} 

= -2tPl p (d u )Pj 2Y (d u )(dZA p , p (n)d^'A Yl (u) + A^ 7 (u)fl2ag'A^(u)) 



2t Sym 



Plp(d u )dZPl 2 y(d U3 ) u lim ui (pJ 7 ,(9 Ul )(A^ a2 (u 13 )^A 7 , 7 (u 12 ) + A^ 7 (u X2 )< Ay aj (u 13 )) 



+ P 7 T 7 , {d U2 ) {A p , a2 (u 23 )d£ A y(s (u 21 ) + A p , p (u 2 i Ay Q2 (u 



13. 



(A.15) 



We have used that d@ A 7 / Q2 (u) is odd, and several times transversality i.e. A a p(u) = P^ j3 (d u )R(u) (this does not 
assume any symmetry) hence Pj iY (d u )Ay a2 (u) = A aiaa (u). 



Appendix B. 1-loop FRG equation in Fourier space 

To one loop one must first solve 



9tCi 3 1 ) a2a3 (k 1 ,k 2 ,k 3 ) = --Sym[P ai ^ 7 (k 1 ) ]T ^^(p.q,^,^)] , (B.l) 

p+q=kj 

where again, here and below Sym[...] means symmetrization w.r.t k" 1 , kj 2 , kj 2 (simultaneous permutations of points 
and indices). Here P a; ^ 7 (k) = ik' 3 P ( J y (k) + ik 7 P ( ^ g (k) for NS and P a; ^ 7 (k) = ik Q <5^ 7 for Burgers, and one uses the 
Gaussian approximation 



C, ^7a 2 a 3 (P' t l' k2 ' k3 ) = ^P+q <5 k 2 +k 3 A ( 3 7 (p)A Q2Q3 (k 2 ) 



+ ^ P +k 2 ^q+k 3 A^ Q2 (p)A 7Q3 (q) + ^p+k 3 <5q+k 2 A^ Q3 (p)A 7a2 (q) 



(B.2) 
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We assume that A Q( g(— k) = A Q!| g(k). The first term vanishes when multiplied by P aii| g 7 (ki). One finds, using 
symmetries 

^^(ki.ka.ka) = « kl+kB+k ,C'^ a3a 3(k lj k 2 ,k3) (B.3) 
^^(ki.ka.ka) = -t(P Ql ^(k 1 )A^ 2 (k 2 )A 7Q3 (k 3 ) 



P a2 ^ 7 (k 2 )A^ 1 (k 1 )A 7Ct3 (k 3 ) + P Q3 ^ 7 (k 3 )A^ 1 (k 1 )A 7Q2 (k 2 ) .(B.4) 



This expression is already symmetric and does not need symmetrization anymore. 
The 1-loop equation is obtained by inserting this result into 



<9 t A QlQ2 (k) = - Sym P Ql(3 <y(k) £ C^, 0a (p,q, -k) 



p+q=k 



(B.5) 



= fSym P Qll3 -y (k)A 7Q2 (k) £ P ( g/ / g 7 (p)A i8 y (q) + P ai(8 'Y (k)A 7C(2 (k) £ Py / g 7 (q)A i8/9 /(p) 

p+q=k p+q=k 

- P Ql( 3- y (k)P ct2(37 (k) A /3/3'(p) A 77'(q) . 

p+q=k 

where here Sym[...] means symmetrization w.r.t. ai,a 2 . Symmetrization finally yields the general 1-loop equation 



(k) =t( (P Ql(9 ,y(k)A 7Q2 (k) +P Q2(3 -y(k)A 7Ql (k)) J2 ^'/3 7 (p)A / 3 7 '(q 



(B.6) 



p+q=k 

-P 0lj9 /y(k)P aaj 9 7 (k) A /3/3'(p) A 77'(q)J > 

p+q=k ' 

from which the 1-loop FRG equation for NS and Burgers can be retrieved. 

For NS one has k Qi A aijQ , 2 (k) = for i = 1,2. One checks on (B.6) that if A is transverse at a given time t, 
it remains so, i.e. the r.h.s. is automatically transverse. For N = 2 this implies that A Q| a(k) = P^g (k)A(k), but 

this is not true for N > 2. The general form is A Q( g(k) = X)fj=i e a(k)e^(k)Aij(k) where the e l (k) span a basis 
orthogonal to k and Ay (k) is a symmetric matrix. 

For simplicity we consider here the subspace A Q( g(k) = P^g (k)A(k). One finds, using Mathematica 

J2 P Ql ^y(k)P Q2 , 7 (k)A^(p)A 77 ,(q) = ^CWJ It^I^jp 

p+q=k q n MJ 



1 



- (k-q)^-O +^(JV-2)fc" ^ + (k-q)" A(q)A(k-q) 



(B.7) 



At this stage this was obtained by: (i) symmetrizing w.r.t. q — > k — q, and (ii) assuming that the result was 
proportional to P^ lQ2 (k), and then contracting with Pj ia2 (k) (or 5 ai _ a2 ). For N = 2 there is no loss of generality, 
and the sum over momenta can be discrete, while for N > 2 this holds only for isotropic turbulence A(k) = A(fc) 
and in the limit of an infinite box where the sums become integrals. 
Next one finds by the same method 



(p Ql(3v (k)A 7Q2 (k) + P Q2(3V (k)A 7Ql (k)) J2 P /3'/j 7 (p)A /3r (q) 



(B. 



p+q=k 



N 



k 2 q 2 - (k • q) 2 



k 2 q 2 (k-q} 2 



(k 2 - q 2 )((k- q ) 2 -q 2 )+(N-2)k 2 (k- q ) 2 A(q)A(k) . 



This yields the 1-loop FRG equation 



uw 2t ^ fcV-(k-q) 2 j pi 
^ A ( k ) = — E fcV(k-q)3 



(k - q) 2 - q' 



(k 2 - q 2 )((k - qf - q 2 ) + (N - 2)k 2 (k - q) 2 A(q)A(k) 



N- 



We note Kraichnan's conventions, 
4/c 2 a k , P ,q = P ai p' 
2fc 2 6 k , p , q = P cjm (k)P, bc (p)pT b (q) . 



(N-2)k 2 [q 2 + (k- q y\ A(q)A(k-q) 



4fc 2 a k , p , q = P Ql/ 3 V (k)P ai/ 3 7 (k)Pj /J ,(p)P 7 r y ,(q) 



(B.9) 

(B.10) 
(B.ll) 
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This corrects a misprint in equation (VII-2-7) of [6]. The various symbols satisfy 



2ak,k-q, q = -4fc 2 a kik - 



[k 2 q 2 - (k • q) 2 ] r 



&k,k-q : q = — 2/j 2 6k,k-q 



»k,k- 



q - q_ k 2 q 2 (k-q) 2 
[k 2 q 2 - (k ■ q) 2 ] 
q ~ k 2 q 2 (k~q) 2 

&k.k-q.q + &k,q,k-q • 



{k-qf-q 2 ) +(N-2)k 2 (q 2 + (k- q ) 



(k 2 -q 2 ) [(k-q) 2 -q 2 +(JV-2)^(k-q) 



The FRG equation can thus be written in various forms, 

It 

-j 5k,k- q , q A(q)A(k - q) - 6 k .k-q. q A(q)A(k) , 



9 t A(k) 



N — 



(B.12) 
(B.13) 
(B.14) 

(B.15) 



as well as the form given in the text. 

Note that &k,k- q , q i s n °t invariant under q — > k — q; cik.k.o = flk,o,k = ^k.k.o while &k,o,k = 0. These properties 
imply that the coefficient B in the expansions in the text is zero. Further k 2 q 2 — (k-q) 2 = k 2 p 2 — (k-p) 2 = p 2 q 2 — (p-q) 2 
if k = p + q, hence this term is already symmetric under q — > k q. 



Appendix C. FRG equation for N = 3 periodic flows 

A convenient parameterization of a general TV = 3 divergence-less velocity correlation matrix, i.e. such that 
k a A a p(k) — 0, (with mirror symmetry) is 

A Q/3 (k) = ^ Eai^k-yEpiskg A t (k) 



^>,(k)(l-k 2 )-A Q (k)^ -k a k? ^A^-A^-A^j 



(C.l) 



In coordinates this is 



A xx (k) = k 2 A z (k) + k 2 A y (k) , A xy (k) = A yx (k) = -k^A^k) 



(C.2) 

and similar for circular permutations. The semi-isotropic case A Q ^(k) = (k)A(k) corresponds to Aj(k) = A(k), 
and is fully isotropic when Aj(k) = A(|k|). For a periodic flow with a cubic lattice symmetry we expect that 

A x (k) = A(k !/ ,k z ;k x ) , A„(k) - A(k x ,k z ;k y ) , A 2 (k) = A(k x , k v ; k 2 ) (C.3) 

where A(ki, k2; k$) is a symmetric function of its first two arguments. 

We have derived FRG equations for the Aj. They are of the form dA z = J2i,j= x y z 2? Aj(q)Aj(k— cfifij (k, q) 
where the fij(k, q) are quite complicated functions of k x ,k v ,k z ,q x ,q v ,q z \ we have not tried to solve them. 



Appendix D. Generating functional approach and diagrammatics 

In this appendix we explain how the FRG equations can be derived for the decaying Burgers, NS and SQG equations 
within the Martin-Siggia-Rose formalism, using a small-time expansion. We start with the Burgers equation, with 
an initial condition v u q = w u at t = + . This is equivalent to 



d t Vut + 2 v u v ut = ^V 2 v ut + 6(t)w u , 



(D.l) 



with v ut = for t < 0, i.e. a forcing which acts only at time zero. We then introduce the generating functional 
e — S[v u t,Vut] f or the velocity correlators, in the usual way, which leads to the dynamic action 



S 



v u t<9tv ut - ^v ut V 2 v ut - -(V„ ■ Vut)(vui) J 



v u qw u 



(D.2) 



v u t is the response field, and the path integral should be evaluated with v ut=0 - = at the boundary. Since 
normalization of the path integral is one, the generating function for averages over the initial conditions can be 
computed from the dynamical path integral with action: 



S' = 



u,t>0 



v u t<9 t v ut - ^v ut V 2 v ut - -(V„ • v ut )(v uf ) 



vS A° (u-u')v&o. (D.3) 



Functional renormalization- group approach to decaying turbulence 



24 



Di 





Figure Dl. The diagrams given in Eqs. (D.7) to (D.9). 



where (w^w^,) ~ A° /3 (u — u'). We use the following graphical representation. The vertex corresponding to the 
cubic nonlincarity is depicted by 

P 



The response function in the limit v — > reads 

a B 

e(* 2 - ti)* a/3 = » . 

1 1 t 2 

The dashed line denotes the 2-point velocity correlator at t = 0, 



(D.4) 



(D.5) 



(D.6) 



We now switch to the Fourier representation. The 2-point velocity correlator written in Fourier space A^g(k) = 

( v kt v -kt) t° 1-loop order is given by the diagrams given on figure Dl. The corresponding expressions with 
combinatorial factors arc 



Di = A^(k), 

D 2 = ^k Q k^A0 5 (q)A° <5 (k + q), 
q 

A3 = -t 2 k Q k /3 ^A° 5 (q)A° 5 (k). 



In real space, the sum of the diagrams can be written as 

A^(u) = A^(u) - ^d a d 3 [A° 5 (u) - A°,(0)] 2 . 

To compute the /3-function we take the derivative with respect to t 

9 t A^(u) = - td a df,[A° 7S (u) - A° d -(0)] 2 , 

and substitute to one loop Aj^g = A^g. As a result we obtain the FRG equation (18). 

We now generalize the method developed above to the Navier-Stokes equation. It is convenient to put all 
derivatives in the cubic vertex on the response field and write the dynamical action as 

v u t dt v ut - v v ut V 2 v u4 - A vi t vZ t (PL (d u ) v£ t J 



(D.7) 
(D.8) 

(D.9) 

(D.10) 
(D.ll) 



S 



vZoKb(u-u')v p u , , (D.12) 



where the last term imposes the initial condition similar to that in action (D.3). The nonlinear cubic vertex can 
then be written as 

IB' 



fp' 

a a — * 



(D.13) 



The two-point function Alg(k) = (v£: t v_ kt ) to second order in A°(k) is given by the diagrams in figure D2. The 
tree-level diagram Si reads 

Si = A° a/3 (k). (D.14) 
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There are six 1PI one-loop diagrams which can be split into two groups. The first group gives 

S 2 + S 3 =t 2 P^(-k)Pl(k)k s k p Y, [A« CT (q)A« p (q + k) + A^(q)AL(q + k)]. (D.15) 

q 

Expanding the projection operators we find 

P^(k)[S 2 + S 3 ] =t 2 J2 a k ,k- q , q A (q)A°(k + q), 
q 

where we used that A° (3 (q) = Pj /3 (q)A°(q) and a is given by equation (B.12). Analogously we obtain 

S 4 + S 5 + S e + S 7 = -t 2 I aT (k)A%(k), (D.16) 
where we have introduced 

i a Ak) = [Pl- t (k)P^ T (k + < l )pjAq)M<i)^(K + ^) 

q 

+ P Q T 7 (k)P 7 T (k + q )Pj CT (q)A( q )k CT (k r + q T ) 
+ P^(k)P 7 T T (k + q)Pj cr (q)A(q)k r (k„ + q M ) 

+ P£(k)J>£(k + q)PX(q)A(q)k 7 (k r + q T )] . (D.17) 
Expanding the projection operators one sees that 

P^(k)[5 4 + S 5 + S 6 + S 7 ] = -t 2 6 k .k- q . q A (q)A°(k), 

q 

where b is defined in equation (B.13). Taking the derivative with respect to t we obtain the FRG equation (26). 

The FRG equation for the SQG equation can be obtained in the same way. The dynamic action corresponding 
to equation (101) with the initial condition imposed at t = is given by 

S = [ {f ut d t T ut - isf ut V 2 T ut + f ut z ■ [VVw x VT ut ]} - I [ T u0 A° (u - u')f u , . (D.18) 

Ju,t>0 1 > I J u>u i 

The SQG vertex can be written in Fourier space as 

f f ut z • [Wu« x VT u4 ] = J f ut e afl dZK u - u ,T uH dPT ut 

J u,f J u,u' .t 

= - I K n _ u ,T u , t e a pdZ[f ut dPT ut ] = - /"^f_ k _ qt e Q(3 q Q k (3 ( 7 - a r qt T kt , (D.19) 

Ju,u',t Jt k q 

where we have introduced the inverse Fourier transform 

K u := FT-^ q ^ • (D.20) 

The 2-point function Aj.(k) = (T k tT_kt) to second order in Aj, is given by the diagrams in figure D2. The 
corresponding expressions read 

Si = A°(k), (D.21) 

^ + ^3 = ^ 2 E A T(q) A T(q + k)[qxk] 2 { g - 2a -2 (Z - a |q + kr a + |q + kr 2a }, (D.22) 
q 

S± + S 5 + S 6 + S 7 = ~t 2 J2 A T(q)A^(k)[q x k] 2 {q- 2a - q- a k- a - q- a \q + k|- Q + £T Q |q + k|- a } . (D.23) 

q 

Taking the derivative with respect to t and reexpressing the bare disorderA^ in terms of the renormalized one A^ 
we obtain the FRG equation (104). 
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Figure D2. Diagrams contributing to the FRG for the NS and SQG equations. 



Appendix E. Distance geometry for FRG equation 

Appendix E.l. Navier- Stokes 

We derive the measure used in [6], e.g. equation (VII-2-9). Following [108], appendix A, the integral of a function 
/(xi, . . . ,x n ) which depends only on Uij := Xj ■ Xj can be written as 



Hd N Xi f( Uij ) 



Sn 



'JV-n+1 



/([««]), 



(E.l) 



where the domain of integration is such that the scalar-products can be realized in iV-dimensional space, and 
5tv = 2n N ^ 2 /T(N/2) is the area of the unit sphere. Here we need n = 2 with k = x 1; p = x 2 , k + p + q = 0, 
k = |k|, q = |q|, p = |p|. One finds, for arbitrary N, 

■,2\H^ SnSn-i ,,. 2 2 



| d^fed^/K) = J dfc 2 dp 2 d(kp)(fcV - (kp) 2 )^^^- 



SnSn-i I dkdpdqkpq 

I A 



/(fc 2 ,?>W), 



(E.2) 



(p+<7~ fc) (k+p—q) (k+q—p) {k+p+q) _ 

where the triangle symbolizes the realization of the triangle inequality. We have used that k • p = i(p 2 + k 2 — q 2 ) 
hence dfc 2 dp 2 d(kp) = idfc 2 dp 2 dg 2 , as well as, 



k 2 p 2 - (kp) 2 = -{k +p + q){p + q - k){k + p - q)(k - p + q). 



(E.3) 



The domain of integration is plotted in figure El. A non-trivial check is to suppose that / is independent of p, and 
to do the p integration. With the domain in figure El, the two cases q < k and q > k have to be distinguished. 
For N = 3, both (!) give J pdp = 2kq, which result in two independent 3-dimensional integrals (with correct 
normalization) for q and k. 

Since this is valid for any function of k, we may rewrite (26) for the isotropic case and any N, replacing 
E 9 f q = I jgSr, using SA(K) = J d N k^=^6A(k), which yields 

r>_y Q c OG c OO o 

LI OjSf—1 I I - ~ * -Ofc 



d t A(k) 



N - 1 (2tt)" Jo 



dq / dpQ A (k,p,q) 



k N pq 



i(fc + p + q){p + q- k){k +p- q)(k -p + q) 
e A (k,p,q) = Q(k+p> q)e(k + q>p)0(p + q> k). 



[A(g)A(p) - A(g)A(fc)] , 



[(k 2 -q 2 )(p 2 ~ q 2 ) + (N-2)k 2 p 2 ], 



(E.4) 



(E.5) 




The domain of integration is plotted in figure El. The distance geometry can be parameterized by p = f (s + t) 
|(s - t) with 



and q 



dp dqQ A (k,p,q)f(k,p,q) 



fc2 
2 



ds 



dtf(k,-(s + t),-(s-t)). 



(E.6) 



ft is immediate from (E.4) that J °° dfc/c w ^ 1 A(fc) is conserved, using the symmetry under the exchange p f-> fc. This 
implies energy conservation for all AT. For N — 2 one checks that enstrophy is conserved, i.e. J °° dfc fc 3 A(fc) can be 
brought to the form ~ J k Bk pq (^ — j-)A(p)A(q) which leads to a factor of (fc 2 — p 2 )(k 2 — q 2 )(p 2 — q 2 ) times a 
symmetric function of p, q; nence it vanishes. 

Violation of energy conservation necessitates a divergence in the integrals for large momenta so that the 
operations involved in the symmetrization, e.g. change of order in integration, are no more valid. For N = 3 one 
sees that at £2 = 1, for fixed s and t there is a logarithmic divergence J dk/k at large fc, none for £2 > 1, and a 
relevant one for C2 < 1- The momentum space- integrals are therefore no longer well-defined for C2 < 1- 



Appendix E.2. Burgers 

One finds a similar expression for Burgers 

SN-lt 



d t A(k) 



4(2tt) 



N 



2-N 



dq / dpQ A (k,p,q)- 



rfc 2 



pq 



(p+q—k) (k+p—q) (k+q—p) (k+p+q) 



(p 2 + q 2 - k 2 ) 2 A(p)A(q) - p 2 (k 2 + q 2 - p 2 ) 2 A(q)A(k) 



(E.7) 



Appendix F. Short-distance expansion of A(u) 
and surface quasi-geostrophic turbulence 



Amplitudes B and C for Burgers, Navier Stokes 



In this appendix, we calculate the necessary integrals for the short-distance expansion of A{u) for Burgers, Appendix 
F.l, Navicr-Stokcs, Appendix F.2 and quasi-geostrophic, Appendix F.3. For simplicity of notations, we use 

b = N + ( 2 . (F.l) 



Appendix F.l. Burgers 

To get the coefficient C in (58) we need to compute 

[q- 



m = h 2 [ [ ^-^ G{ q)G { k - g) = h 2 [ k ( ^-f )2 



G 



2 +q 



G 



2 q 



(F.2) 



and set G(k) = l/k b . We note gencrically / an integral where such replacement is performed, while J denotes 



the same integral with IR cutoffs. Indeed, the second integral in (58) behaves as 



,N-bu2 



k 2 /{k 2 



ra 2 ) b / 2 and 



contributes only to B and D (equivalent to the statement that f q q a = in dimensional regularization) . While 
the calculation of B and D (for each integral) depends a priori on the IR details of G(fc), the coefficient C can 
be obtained by the method of analytical continuation on the first integral only. This integral is both UV and IR 
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convergent upon inserting G(k) = l/k b for N/2 < b < N. Its expression is then continued for b > N to get C. The 
cancellation of B between the two terms is easy to show. One has 

(^-<7 2 ) 2 k 2 



k 2 

m - * T 



1 



q i 2 



| + q| b+2 l|-q| 6+2 2r(6/2 + i) 5 



U>0 



t b/2 t b/2Q 2 \ / e -(t 1 +t 2 )(if+ g 2 )+ V (^- 3 2 )-(4 1 -t 2 )k., 

1 2 v | u= o / 



(47T) 



-AT/2 



2 r(6/2 + l) 2 



ii /2 4 /2 9 2 L =0 (ii + * 2 + v )-^ e ^+W^ 



fi>0 



AT/2, 



-*«1«2 fc 2 



(47T) 



-N/2 



U 2 

2 4r(6/2 + l)U I >o 



,6/2,6/2/, , , \-4-^ -Pxi-k 2 

t{ t 2 ' (ti +t 2 ) 2 e *i+ t2 



iV 2 (ti + t 2 ) 2 + 27V (ti + i 2 ) (*2 + ti(l - 2k 2 t 2 )) + 2fc 2 <J - t x t% + t\- t\t 2 {\ - 2fc 2 t 2 ) 
Introducing s = t± + t 2 , t\ = su, t 2 = s(l — u) with d<idt2 = sduds one gets 



(F.3) 



V ; 2 4r(6/2 + l) 2 y 



du 



,-1+6-f 



i5>0 



- u)) b/2 \ \N 2 + N(2 - 4fc 2 s(l - u)u) 



2k 2 s 



l-2(l-u)u(2-fc 2 s(l 
The integration over s can be performed, if b > N/2, leading to 

m . - ^Q'^^' ^-jf' d. Mi - .)]«*=» [» - 4*1 - »)« 

- 4iV(l + JV)(1 - u)u - 26(1 - 2(1 + 2N)(1 - u)u) 
This integral converges only for b < N . where it is 

_ ^,2+^-25 _ -2+6-2AT_i^ _ 2b ( N _ 



3 — su(l— u)k 



(F.4) 



7(/c) = 



a+6 _ 2 v [6 2 - 26(7V - 1) + iV(JV - l)]r(6 - f )r(l (JV - h)) fc2+jV _ 
Y{b/2 + lfT{\{l + N -b)) 



(F.5) 



(F.6) 



This identifies C = C(b, N) for b < N. For b > N the correct calculation requires regularization by a mass m 2 , and 
leads to J(k) = B\k 2 ~ b + Ck 2+N ~ 2b + Dk~ b . While B\ is cancelled by the second integral in (58), the expression 
of C remains equal to the analytical continuation of (F.6). 



Appendix F.2. Navier-Stokes 

The non-linear term in the FRG equation is the sum of two integrals, 
k 2 q 2 - (k • q) 2 



(k 2 - q 2 )[(k - q y-q 2 \+(N~ 2)fc 2 (k - q) 2 G(q) G(k - q) - G(k) 



J(fc) (N-l)Jq F g 2 (k-q)2 

If we replace G(p) by p~ b , we see that the first integral is both UV and IR convergent for N/2 < b < N, while the 
second is nowhere convergent. It is either UV divergent (for b < N) or IR divergent (for b > N) at q w (but not 
at q « k). It is thus convenient to split J(k) into two parts: 

J(k) =J(k) + J c (k), (F.7) 

fc 2 g 2 - (k ■ q) 2 



Jc(k) 



J(k) 



-G(q) G(|k-q|)-G(*) 



k 2 q 2 — (k • q) 



2 . 



q 2 - k 2 - (k - q) 2 G(q) G(|k - q|) - G(k) 



(F.8) 
(F.9) 



(7V-l)7 q P(k-q)2 

It is easy to see that the second integral in J c (k) contributes only to B and D but not to C, while I(k), the same 
integral as J(k), replacing G(p) by p~ b , is now both UV and IR convergent for N < b < N + 2. One thus has 



I(k) 



2k 2 (N-l)J tut2>0 



t 

1 i 



6/2-1^/2 



t^k-^ 



r(i + |)r(|) r(|) 



B(h,t 2 ) = [k A + 2k 2 (d tl +d t2 ) + (d tl -d t2 ) 2 ] (-k 2 + d t2 -d tl ) , 



using that k q — (k ■ q) 



fc 4 - 2k 2 (q 2 + (k - q) 2 ) + (q 2 - (k - q) 2 )' 



Using 



-tl<J 2 -t 2 (k-q 



) 2 = (47T)- N / 2 (t 1+ t 2 )- N / 2 e- k2 ^ 



(F.10) 
(F.ll) 

(F.12) 
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(F.13) 



Introducing s = ii + iaj £1 = su, t 2 = s(l — u) with dtidt2 = sduds one gets 



I(k) = 



2k 2 



(47r)^/2r(i + |)r(|) Jo 



du 



>o 



S-l+f-f [ U (l _ U )]l |l _ s -6/2(l _ u )-6/2r^l + -) k- 



— k su(l—u) 



The integral over 5 can be performed in the first integral for b > N/2 and for b > N in the second. This leads to 



i{k) = 



(4^/ 2 r(i + !)r(|) 



b--)t 



.2-2Z.+JV 



du [u(l — u)]~ 



r(i + |)r(^) 
r(6-f) 



Both integrals are convergent for b < N + 2 and one finds I(k) = (Ci + C^fc 2 2b+JV with 
2 b - jV 0Fbr(6- f )r(|(2-6 + iV)) 



Co 



2(4^) w / 2 r(i + |) 2 r(i(3 - b + n)) 



{4ir) N / 2 sin (i(6 - iV)7r)r(^)r(2 - | + iV) 
We now compute the integral associated to J c (k) 



(F.14) 

(F.15) 
(F.16) 
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The integral over s can be done for b > N/2, and the one over u for b < N. It gives 

2 b ~ N (N - l)y/nT(b - f )r(±(iV - 6)) 



Ic{k>) — C c k 



2-26+ JV 



a = 



(4^/ 2 r(i + |)r(|)r(i(i - 6 + n)) 



The total contribution is thus 

J(jfe) = Bfc 2 - fc + Cfc 2 - 2fc+w + ^fc- fc 
C - Ci + C 2 + C c , 



(F.17) 



(F.18) 
(F.19) 



where G\ + C c is the total contribution to C from the first integral, and Ci the contribution of the second integral. 
Of course one also shows B = 0, as a result of the cancelations. 

Appendix F.3. Surface quasi- geostrophic turbulence 

Let us study the nonlinear term in the FRG equation for the SQG turbulence given by 

d 2 q 



J(k) 



C 7^ 

J (2^ [q X k]2[q ~ a -P~ a )&T(q)[(q~ a - P - a )A T (p) - 2(q- a - k- a )A T (k)], (F.20) 



where p = k q. We now assume for the correlator Ar(fc) the form 
It is related to the velocity correlator A(k) by 



/3 = 4-2a + C 2 , 



A(fc) 



fc2+C2 



(F.21) 
(F.22) 
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Figure Fl. Blue solid curves: Locations in the (a, /3) plane, where F(a,/3) given in equation (F.27) vanishes. The 
two straight lines are given by fi\ = 1 + ^(2 — a) and /32 = 4 — a. Orange dashed line: /3 = 6 — 2a. 



Keeping in mind that the case a = 2 corresponds to the NS equation in N = 2, we expect that f3 is close to 2 near 
a = 2. We thus want to compute 

d 2 q 
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(2n)' 
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Using that 



we obtain 
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where B{t 1 ,t 2 ) = fc 4 + 2fc 2 (<9 tl + 9 t J + (<9 tl - dt 2 f . Using equation (F.12) for N = 2 we find 
I{k) =F(a,/3)fc~ 2(Q+/3 - 3) , 

£^ r _ £ _i_ o\ y^2 a+l3 T (-f 



(F.26) 
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(F.27) 



One checks that for a = 2 one recovers C(N = 2,(2) = F(a = 2,(3 = C2) where C(N,( 2 ) is given by equation 
(61). As we already know this means that the limits (3 — > 2 and a — > 2 are not exchangeable without an IR cutoff. 
Indeed, we have lima^ F(a = 2,(3) = — 1/(4-71") and not — l/(87r), which we expect in the presence of an IR cutoff. 
However, we find that if one keeps o ^ 2 infinitcsimally close to 2, and then takes the limit (3 -+ 2, one does find 
— l/(87r). This means that one needs to keep a > 2, and that lim a _ >2 F(a, (3 = 2) = — l/(87r) is the correct limit. 
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This result leaves two options for the large-fc behavior of possible fixed points parameterized by /?: 

(i) /3 = 6 — 2a leads to the nonlinear term in the FRG equation which for large fc asymptotics is also ~ Fk~@. 
Thus, it can be balanced by the rescaling terms and may have a self consistent solution even if F is not zero. The 
tail is then given by Ax(k) ~ l/fc 6_2a . In this case the tail A(fc) ~ 1/fc 4 for the velocity remains independent of a. 

(ii) there is a value of f3 < 6 — 2a where F(a,f3) vanishes. All possible values are shown on figure Fl for 
— 2 < a < 2. For instance we find that F(a, (3) vanishes for /?i = 1 + |(2 — a) and fa = 4 — a, which cross at a = 0. 
However, there arc other values. 

If we ask that the function £2(0) for the physical fixed point is continuous, it follows from our analysis for the 
2D NS equation (where £2(0 = 2) = 2) that likely values are (2(a) — 2 or £2(0) = a for a > 0. However more work 
is clearly called for, in view of these results, to study possible fixed points as a function of a. 



Appendix G. Small-fc expansion for Navier-Stokes in dimension N = 2 



In this appendix, we calculate the small-fc expansion of the nonlinear term in the FRG equation for the 2-dimensional 
decaying Navier-Stokes turbulence. This will allow us to find in a self-consistent way the small-fc behavior of the 
FP solution A*(fc). To this aim, consider the non-linear contribution to the flow-equation in the distance geometry 
representation given by equation (71). It is useful to symmetrize it in t — > —t: 
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vV-l)(l-* 2 ) 
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{[((* - t? 4) £(§(«-<)) - {(s + t) 2 -4) A(f (s+i))] A(fc) + 4stA(l(s-t))A(l(s+t))} 



We now want to expand equation (G.l) in small k for an arbitrary function A(fc). In Sec. 6.2.2 we already discussed 
that the expected behavior of the fixed-point solution A*(k) at small k is A(fc) ~ fc™ -1 with n = 3. Let us for 
the moment consider a more general class of functions with a finite A(fc = 0) and A'(fc = 0). It is not possible to 
expand the integrand of (G.l) in small k, since this gives integrals diverging at large s. Instead, one can rescale s, 
by defining s = 2q/k, and only then expand in k. This allows one to integrate term by term over t € [0, 1], 
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We have used integrations by parts and the large-A; behavior (73) which suggests that lim/ ; _j. 00 k 2 A(k) = 0. Note 
that the expansion of A(q) in the integrands of the first line of equation (G.2) can not produce terms of order fc 4 . 
Expanding (G.l) further to order fc 4 , we find 
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Together with (G.2), this yields 
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(G.4) 



Note that at least the first few terms of the expansion (G.4) depend mainly on the integral properties of A(fc) and 
not on its small- fc expansion. 



Appendix H. Asymptotic large-fc behavior for Navier-Stokes in N = 2 



In this appendix, we calculate the asymptotic large-fc behavior of the non-linear term in the flow-equation for 2D 
Navier Stokes. We start from the rescaled dimensionless version of (26), setting 

A « - why- w 

Thus we need to compute for N = 2 the convergent integral 
d 2 q fc 2 q 2 - (k • q^ 2 - 
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(H.2) 
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After rcscaling — >• qk and m — > mik it becomes 
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Changing variables to t in such a way that cos(</>) 



t 2 + 



Y and sin(0) — > -prpj with Jacobian t ^^_ 1 we arrive at 
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The integral over t has to be taken independently for < q < 1 and 1 < q < oo so that 6A(k) = 5iA(k) + 52A(k). 
Introducing A = y^rof + (g — l) 2 and B = \J to 2 + (q + l) 2 one can write these integrals in the following form 
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Integration over q and combining both terms gives 
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with 
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(H.9) 



Expanding in small 21, i.e. large fc/m, we obtain 

/(*) = 1 + a 2 f 8 log(x) - - x 4 (32 logOr) + f ) + O (x 5 ) . 

Note that only the leading term is universal, while the higher ones depend on the regularization by to introduced 
in equation (H.l). This implies the leading-order term — given in equation (80). 



(H.10) 




Figure II. Left: Check for the precision of the fit for f(t) defined in Eq. (1.2). Right: The guessed fixed-point 
A gU ess(£:) (blue, thin line), and the numerical solution A(fc) defined in Eq. (84) (red, fat line). See figure 2 for a 
log-log plot of A(fc). 




Figure 12. Top left: <5A(fc) (green solid line with superimposed blue dots), 2A(fc) + fcA'(fc) (solid red line) and 
diA(k) (dashed black line). Top right: The same 3 terms, multiplied by (f + fc 2 ) 2 . Bottom: plot of the ratio 

— . — ^ A ( fc ) which must be 1 at the fixed point. 

2A(fc)+fcA'(fc) r 



Appendix I. Numerical solution for the fixed point in dimension N = 2 

We wish to integrate equations (65) and (71) numerically, using ( = 2. If the fixed point is attractive, A(fc) will 
converge against it. Numerically, the problem is hard for several reasons: 

(i) convergence of the integral (71) for large k is slow and imprecise. 

(ii) one needs high precision for A(fc), which respects the conditions (74) and (78) at small k, and the asymptotic 
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form (80) for large k. 

(iii) the r.h.s. of equation (65) must be calculated numerically at many different points; we will use 108 points (see 
table II). 

(iv) one easily runs into numerical instabilities, when using a spline-interpolation through all points, or a polynomial 
fit of degree as 20, which is necessary to represent faithfully the data-points. 

In order to circumvent these problems, we do the following: Write the ansatz (84)-(85). The guess in equation (85) 
was obtained by (i) imposing the correct asymptotic forms, (it) trying to optimize consistency relations as (G.4). 
Since we have no small-fc expansion with rapidly decaying coefficients, the final criterion (iii) for (85) was a r.h.s. 
in the flow-equation which is small in the intermediate regime of finite k. 

Our information is stored in A cor (fc), spaced as the first column in table II. It is updated via 

A cot (k) -> A cor (fc) + Ktd t K{k) , (1.1) 

at the given values of k. The quantity d t A(k) is given by equations (65) and (71). Note that while we use the 
distance geometry formula (71) for the correction since it is numerically more efficient, we have checked that the 
result is the same when replacing equation (71) by equations (66) and (67). The step-size k is finally reduced to 
k = 0.005. 

To circumvent the numerical problems mentioned under (iv) above when obtaining the rescaling term in (65) 
and the numerical integral in (71), we approximate A cor (fc) by the smooth function 

a_w-^ (of 

m(k) 
u2 , i\6/a 



m (k) := ± t^T— (1-4) 



(k 2 +_l 

k 2 

The function fit) is defined on [0, 1], bounded (of order 100), and converges to for t to or 1, see figure II. 
In order to produce a smooth fit for f(t), we use the best (i.e. least variance) cubic spline with o equidistant 
points and /(0) = /(l) = /'(0) = /'(l) = 0. o is initially chosen to be 13, and then increased up to 20. Using 
splines of this relatively low order (compared with the number of data points) effectively filters out numerical noise. 
(Note: A polynomial fit of the same order is not adequate, but generates a dynamic instability. To further increase 
precision, a spline with variable switch-points could be used, o = 20 is our upper limit due to RAM problems in the 
implementation of the integration routines, which split the integrals in pieces for each spline part.) Before updating 
A cor (A:) via (1.1), the table of stored values for A cor (fc) is replaced by the approximation (1.2). This is necessary for 
consistency. On the left of figure II, we show our final result of the function f(t), and the points from which it is 
constructed. One sees that no numerical artefact is present. 

Our final result for A(k) is shown in a log-log plot on figure 2. A linear plot is presented on the right of figure 
II (fat red line), as well as our initial guess A guess (fc). Note that since there exists a redundant mode (the choice 
of fe m ax at which E(k) = fcA(fc) is maximal, one could get these two curves closer. Numerical values for A(fc) are 
given in table II. The precision can be inferred from figure 12. It should be few percent (relative precision) for 
k < 5, but the precision decreases for larger values of k. Note however that we know the exact asymptotic form 
A(fc) = 167r/fc 4 . 



Appendix J. Navier-Stokes equation in the limit of large N 



In this appendix we consider the FRG flow for the Navier-Stokes equation (32) in the limit of large N. We start 
from the flow equation in real space (32) to derive an equation for V„A(u): 



tdtVl A(u) = (2 - C) V 2 A(u) + ^ud u [v£A(u)] + V^Al(u) + V*<fA NL (u). 



(3.1) 



The Laplacian is V„ = Nd/dy + 2yd 2 /dy 2 ; rewriting A(u) in terms of r(y), given by equation (31), we obtain 



ViA(u) 



N (2 + N) r"(y) + Ay ((2 + N) r"'(y) + y r (4) (y)) 



(J.2) 



The Laplacian of V 2 1 <SAl(m) is local and has been expressed in terms of r(y) in equation (33). The nonlocal part 
<5Anl(w) is given by equation (34) with the kernel (35). Taking into account that for a large- N expansion in real 
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k 


A(fc) 


k 


A(k) 


k 


A(k) 


k 


A(k) 








0.869074 


30.8943 


2.0025 


18.629 


4.29082 


0.264334 


0.0158116 


0.0131526 


0.9045 


32.5893 


2.05554 


16.0965 


4.44768 


0.222674 


0.0474395 


0.117591 


0.940243 


34.2548 


2.10983 


13.8072 


4.6164 


0.186731 


0.0790817 


0.323878 


0.976321 


35.9175 


2.16544 


11.7798 


4.79873 


0.15576 


0.110748 


0.629261 


1.01275 


37.6048 


2.22246 


10.005 


4.99681 


0.12912 


0.142447 


1.03318 


1.04955 


39.3028 


2.28096 


8.46986 


5.21327 


0.106266 


0.174189 


1.53977 


1.08673 


40.9413 


2.34104 


7.15983 


5.45142 


0.0867232 


0.205984 


2.15575 


1.12433 


42.4538 


2.40281 


6.05898 


5.71548 


0.0700738 


0.237841 


2.8848 


1.16235 


43.782 


2.46637 


5.14416 


6.01089 


0.0559517 


0.26977 


3.72542 


1.20082 


44.876 


2.53185 


4.38304 


6.34489 


0.044037 


0.301781 


4.66717 


1.23977 


45.7066 


2.59937 


3.74686 


6.72728 


0.0340508 


0.333885 


5.68728 


1.2792 


46.2688 


2.66909 


3.21119 


7.17178 


0.0257567 


0.366091 


6.76818 


1.31916 


46.5628 


2.74115 


2.7555 


7.69832 


0.0189527 


0.39841 


7.91167 


1.35967 


46.5915 


2.81574 


2.36475 


8.33709 


0.0134578 


0.430852 


9.12235 


1.40074 


46.3608 


2.89304 


2.02966 


9.13673 


0.00910839 


0.463428 


10.4041 


1.44242 


45.8577 


2.97328 


1.74234 


10.1818 


0.00575792 


0.496149 


11.7583 


1.48473 


45.0426 


3.05668 


1.49586 


11.6362 


0.00328225 


0.529027 


13.1873 


1.5277 


43.8845 


3.14352 


1.28418 


13.8756 


0.00157192 


0.562072 


14.6963 


1.57137 


42.3626 


3.23409 


1.10221 


16 


0.000868729 


0.595297 


16.2907 


1.61578 


40.4659 


3.32872 


0.945681 


18.0514 


0.00052673 


0.628713 


17.9758 


1.66095 


38.2213 


3.42779 


0.810933 


20 


0.000344759 


0.662334 


19.758 


1.70694 


35.7034 


3.53174 


0.694791 


23.5 


0.000177434 


0.69617 


21.6242 


1.75378 


32.9836 


o n A 1 f\ A 

3.64104 


0.59454 


27 


0.000100406 


0.730237 


23.5339 


1.80152 


30.1263 


3.75626 


0.507916 


31.5 


0.0000534954 


0.764546 


25.447 


1.85021 


27.1889 


3.87803 


0.433078 


35.5 


0.0000328852 


0.799113 


27.3277 


1.8999 


24.2391 


4.00711 


0.368428 


41 


0.0000183294 


0.83395 


29.1463 


1.95064 


21.362 


4.14436 


0.312578 


48 


9.68650 xlO" 6 



Table II. Numerical fixed point of equations (65), (71) for A(fc). 



space, both p 2 and k 2 ~ N, we can expand the kernel in 1/N as follows 



(J.3) 



N k 2 + p 2 

The first term in equation (J.3) can be written as a derivative, thus becomes also local 

V 2 u 6A$l(u) = 2(N-l) r'(0) [(4 + IV) ((2 + TV) (n r'" (y) + 6yr^(v)) + ^y 2 r^{y)) + 8y 3 r^(y) 

Note that apart from the factor of r'(0), it has the structure of the last term in equation (33). Let us now evaluate 
the second term in equation (J.3) to leading order at large N, 



(J.4) 



4 
N 



-iku 7„6 



k b R(k)- 



1 



-p 4 R( P ) 



k . p 



k 2 + p 2 
~" I fc 6 jj(fc)p 4 J R(p)e' 



A:./) 



Since J k e R(k) = R(u) = r(y), we have in the large-A/" limit 



-iku ^—tk 



e~ tk R(k) = r(y + Nt) , 

as can be checked by Taylor expanding and using that — k 2 = V„ = Nd/dy to leading order. Hence 



(J.5) 



(J.6) 



l NL \ —j N 

We now put all pieces together, setting £ = (q + (-i/N 



dtr"'{y + Nt)r"{Nt) = —N 4 



dzr"'(y + z)r"(z). 



(J.7) 



Expanding to order 1/N we obtain equation (94). 



Functional renormalization- group approach to decaying turbulence 



36 



References 

[1] A.N. Kolmogorov, On the energy distribution in the spectrum of a turbulent flow, C.R. Acad. Sci. URSS 30 (1941) 301-305. 
[2] A.N. Kolmogorov, A refinement of previous hypotheses concerning the local structure of turbulence in a viscous incompressible 

fluid at high reynolds number, J. Fluid Mech. 13 (1962) 82-85. 
[3] U. Frisch, Turbulence, Cambridge University Press, 1995. 

[4] A.M. Oboukhov, Some specific featues of atmoshperic turbulence, J. Fluid Mech. 13 (1962) 77—81. 

[5] R.H. Kraichnan, The structure of isotropic turbulence at very high reynolds numbers, Journal of Fluid Mechanics 5 (1959) 
497-543. 

[6] M. Lesieur, Turbulence in Fluids, Kluwer Academic Publishers, Dordrecht/Boston/London, third edition, 1995. 
[7] C.-Y. Mou and P.B. Weichman, Spherical model for turbulence, Phys. Rev. Lett. 70 (1993) 1101-1104. 

[8] C.-Y. Mou and P.B. Weichman, Multicomponent turbulence, the spherical limit, and non-kolmogorov spectra, Phys. Rev. E 52 

(1995) 3738-3796. 

[9] R.H. Kraichnan, Small-scale structure of a scalar field convected by turbulence, Phys. Fluids 11 (1968) 945-953. 
[10] D. Forster, D.R. Nelson and M.J. Stephen, Large- distance and long-time properties of a randomly stirred fluid, Phys. Rev. A 
16 (1977) 732-749. 

[11] C. DcDominicis and P.C. Martin, Energy spectra of certain randomly -stirred fluids, Phys. Rev. A 19 (1979) 419-422. 

[12] J.R. Mayo, Field theory of the inverse cascade in two-dimensional turbulence, Phys. Rev. E 72 (2005) 056316. 

[13] L.T. Adzhemyan, N.V. Antonov, P.B. Gol'din, T.L. Kim and M.V. Kompaniets, Renormalization group in the infinite- dimensional 

turbulence: third-order results, J. Phys. A 41 (2008) 495002. 
[14] L.Ts. Adzhemyan, N.V. Antonov, M.V. Kompaniets and A.N. Vasil'ev, Renormalization- group approach to the stochastic navier- 

stokes equation: Two-loop approximation, Int. J. Mod. Phys. B 17 (2003) 2137-2170. 
[15] LTs. Adzhemyan, NV. Antonov and AN. Vasil'ev, Infrared divergences and the renormalization group in the theory of fully 

developed turbulence, Zhurnal Ekspcrimcntalnoi i Tcorctichcskoi Fiziki 95 (1989) 1272-8. 
[16] A.M. Polyakov, Turbulence without pressure, Phys. Rev. E 52 (1995) 6183-6188, hep-th/9506189. 

[17] S. Boldyrev, T. Linde and A. Polyakov, Velocity and velocity- difference distributions in Burgers turbulence, Phys. Rev. Lett. 

93 (2004) 184503, nlin.CD/0305014. 
[18] A.M. Polyakov, The theory of turbulence in 2- dimensions, Nucl. Phys. B 396 (1993) 367-385, hep-th/9212145. 
[19] W. Pauls, T. Matsumoto, U. Frisch and J. Bee, Nature of complex singularities for the 2D Euler equation, Physica D 219 (2006) 

40-59. 

[20] M. Chertkov, A. Pumir and B.I. Shraiman, Lagrangian tetrad dynamics and the phenomenology of turbulence, chao-dyn/ 
9905027 (1999). 

[21] V.S. L'vov, E. Podivilov, A. Pomyalova, I. Procaccia and D. Vandembroucq, Improved shell model of turbulence, Phys. Rev. E 
58 (1998) 1811-1822. 

[22] L. Ts. Adzhemyan, N.V. Antonov, J. Honkonen and T. L. Kim, Anomalous scaling of a passive scalar advected by the navier-stokes 

velocity field: Two-loop approximation, Phys. Rev. E 71 (2005) 016303. 
[23] M. Chertkov, G. Falkovich, I. Kolokolov and V. Lebedcv, Normal and anomalous scaling of the fourth-order correlation function 

of a randomly advected passive scalar, Phys. Rev. E 52 (1995) 4924-4941. 
[24] K. Gawedzki and M. Vergassola, Phase transition in the passive scalar advection, Physica D 138 (2000) 63—90. 
[25] A.L. Fairhall, O. Gat, V. L'vov and I. Procaccia, Anomalous scaling in a model of passive scalar advection: exact results, Phys. 

Rev. E 53 (1996) 3518-35. 

[26] G.L. Eyink and J. Xin, Self-similar decay in the Kraichnan model of a passive scalar, J. of Stat. Phys. 100 (2000) 679-741. 
[27] K. Gawgdzki and M. Vergassola, Phase transition in the passive scalar advection, cond-mat/9811399 (1998). 
[28] V.S. L'vov, I. Procaccia and A.L. Fairhall, Anomalous scaling in fluid mechanics: the case of the passive scalar, Phys. Rev. E 
50 (1994) 4684-704. 

[29] K.J. Wiese, The passive polymer problem, J. Stat. Phys. 101 (2000) 843-891, chao-dyn/9911005. 

[30] L.T. Adzhemyan and N.V. Antonov, Renormalization group and anomalous scaling in a simple model of passive scalar advection 

in compressible flow, Phys. Rev. E 58 (1998) 7381-96. 
[31] D. Bernard, K. Gawedzki and A. Kupiainen, Anomalous scaling in the n-point functions of a passive scalar, Phys. Rev. E 54 

(1996) 2564-2572.' 

[32] D. Bernard, K. Gawedzki and A. Kupiainen, Slow modes in passive advection, J. of Stat. Phys. 90 (1998) 519-69. 

[33] S. Chen and R.H. Kraichnan, Simulations of a randomly advected passive scalar field, Phys. Fluids 10 (1998) 2867-84. 

[34] M. Chertkov and G. Falkovich, Anomalous scaling of a white advected passive scalar, Phys. Rev. Lett. 76 (1996) 2706. 

[35] U. Frisch, A. Mazzino, A. Noullez and M. Vergassola, Lagrangian method for multiple correlations in passive scalar advection, 

Phys. Fluids 11 (1999) 2178-86. 
[36] K. Gawgdzki and A. Kupiainen, Anomalous scaling of the passive scalar, Phys. Rev. Lett. 75 (1995) 3834. 

[37] R.H. Kraichnan, V. Yakhot and S. Chen, Scaling relations for a randomly advected passive scalar field, Phys. Rev. Lett. 75 
(1995) 240-3. 

[38] G. Falkovich, K. Gawedzki and M. Vergassola, Particles and fields in fluid turbulence, Rev. Mod. Phys. 73 (2001) 913. 

[39] O. Gat, VS. L'vov, E. Podivilov and I. Procaccia, Nonperturbative zero modes in the Kraichnan model for turbulent advection, 

Phys. Rev. E 55 (1997) R3836-9. 
[40] B.I. Shraiman and E.D. Siggia, Symmetry and scaling of turbulent mixing, Phys. Rev. Lett. 77 (1996) 2463-6. 
[41] G.K. Batchelor, Computation of the energy spectrum in homogeneous twodimensional turbulence, Phys. Fluids 12 (1969) 11-233. 
[42] D. Bernard, G. Boffctta, A. Celani and G. Falkovich, Inverse turbulent cascades and conformally invariant curves, Phys. Rev. 

Lett. 98 (2007) 024501. 

[43] M.A.I. Flohr, Two-dimensional turbulence: A novel approach via logarithmic conformal field theory, Nucl. Phys. B 482 (1996) 
567 - 578. 

[44] C. Eling, I. Fouxon and Y. Oz, The incompressible Navier-Stokes equations from black hole membrane dynamics, Phys. Lett. B 
680 (2009) 496-499. 

[45] I. Bredberg, A. Strominger Black Holes as Incompressible Fluids on the Sphere, arXiv:1106.3084. 
[46] I. Bredberg, C. Keeler, V. Lysov, A. Strominger From Navier-Stokes To Einstein arXiv:1101.2451 



Functional renormalization- group approach to decaying turbulence 



.37 



[47] M. Muller, J. Schmalian and L. Fritz, Graphene: A nearly perfect fluid, Phys. Rev. Lett. 103 (2009) 025301. 
[48] J. Bee and K. Khanin, Burgers turbulence, Phys. Rep. 447 (2007) 1-66. 

[49] P. Le Doussal, Chaos and residual correlations in pinned disordered systems, Phys. Rev. Lett. 96 (2006), cond-mat/0505679. 
[50] U. Frisch and J. Bee, Burgulence, in New Trends in Turbulence, Springer EDP-Sciences, 2001. 

[51] J. P. Bouchaud, M. Mezard and G. Parisi, Scaling and intermittency in Burgers turbulence, Phys. Rev. E 52 (1995) 3656-3674, 
cond-mat/9503144. 

[52] M. Mezard, Disordered systems and Burger's turbulence, J. Phys. IV (France) 8 (1997) 27-38. 

[53] P. Le Doussal, M. Muller and K.J. Wiese, Avalanches in mean-field models and the Barkhausen noise in spin-glasses, EPL 91 

(2010) 57004, arXiv: 1007.2069. 

[54] P. Le Doussal, M. Muller and K.J. Wiese, Cusps and shocks in the renormalized potential of glassy random manifolds: How 
functional renormalization group and replica symmetry breaking fit together, Phys. Rev. B 77 (2007) 064203, arXiv:0711.3929. 
[55] J.D. Fournier, U. Frisch and H.A. Rose, Infinite- dimensional turbulence, J. Phys. A 11 (1978) 187-198. 
[56] M. Hairer, Solving the KPZ equation, arXiv:1109.6811 (2011). 

[57] D.S. Fisher, Interface fluctuations in disordered systems: 5 — e expansion, Phys. Rev. Lett. 56 (1986) 1964-97. 

[58] P. Le Doussal, K.J. Wiese and P. Chauve, Functional renormalization group and the field theory of disordered elastic systems, 

Phys. Rev. E 69 (2004) 026112, cond-mat/0304614. 
[59] K.J. Wiese and P. Le Doussal, Functional renormalization for disordered systems: Basic recipes and gourmet dishes, Markov 

Processes Relat. Fields 13 (2007) 777-818, cond-mat/0611346. 
[60] P. Le Doussal, Functional renormalisation for disordered systems, in Euro-Summer School on Condensed Matter Theory 2004, 

lecture notes at http://www.lancs.ac.uk/users/esqn/windsor04/program.htm. 
[61] P. Le Doussal, Finite temperature Functional RG, droplets and decaying Burgers turbulence, Europhys. Lett. 76 (2006) 457-463, 

cond-mat /0605490 . 

[62] P. Le Doussal, Exact results and open questions in first principle functional RG, Annals of Physics 325 (2009) 49-150, 
arXiv:0809.1192. 

[63] L. Balents, J. P. Bouchaud and M. Mezard, The large scale energy landscape of randomly pinned objects, J. Phys. I (France) 6 

(1996) 1007-20. 

[64] P. Le Doussal and K.J. Wiese, Size distributions of shocks and static avalanches from the functional renormalization group, 

Phys. Rev. E 79 (2009) 051106, arXiv:0812.1893. 
[65] P. Le Doussal and K.J. Wiese, First-principle derivation of static avalanche- size distribution, Phys. Rev. E 85 (2011) 061102, 

arXiv:1111.3172. 

[66] P. Le Doussal, A. Rosso and K.J. Wiese, Shock statistics in higher- dimensional Burgers turbulence, EPL 96 (2011) 14005, 
arXiv:1104.5048. 

[67] V.S. L'vov, R.A. Pasmanter, A. Pomyalov and I. Procaccia, Strong universality in forced and decaying turbulence in a shell 

model, Phys. Rev. E 67 (2003) 066310. 
[68] P. Le Doussal, What can disordered systems tell us about turbulence, in The Nature of Turbulence, KITP program 2011, 

http:/ /online. kitp.ucsb.edu/online/turbulencell/ledoussal/ 
[69] J. Leray, Etude de diverses equations integrales nonlineaires et de quelques problemes que pose I'hydrodynamique, J. Math. Pure 

Appl. 12 (1934) 1-82. 

[70] J. Leray, On the movement of a viscous fluid to fill the space, Acta Mathcmatica 63 (1934) 193-248. 

[71] J. Duchon and R. Robert, Inertial energy dissipation for weak solutions of incompressible Euler and Navier-Stokes equations, 

Nonlinearity 13 (2000) 249-255. 
[72] K. Gawcdzki, Easy turbulence, arXiv:chao-dyn/9907024 (1999). 

[73] S. Kida, Asymptotic properties of burgers turbulence, J. Fluid Mech. 93 (1979) 337-377. 

[74] J. -P. Bouchaud and M. Mezard, Universality classes for extreme-value statistics, J. Phys. A 30 (1997) 7997—8015, arXivxond- 
mat/9707047. 

[75] A. Rosso, Y.V. Fyodorov, P. Le Doussal, Freezing transition in decaying burgers turbulence and random matrix dualities, 

Europhys. Lett 90 (2010) 60004, arXiv:1004.5025. 
[76] S.N. Gurbatov, S.I. Simdyankin, E. Aurell, U. Frisch and G. Toth, On the decay of burgers turbulence, J. Fluid Mech. 344 

(1997) 339-374. 

[77] L. Balents and D.S. Fisher, Large-N expansion of 4 — e-dimensional oriented manifolds in random media, Phys. Rev. B 48 
(1993) 5949-5963. 

[78] J.M. Burgers, The non-linear diffusion equation, Dordrecht, 1974. 

[79] T. de Karman and L. Howarth, On the statistical theory of isotropic turbulence, Proc. R. Soc. Lond. A 164 (1938) 192-215. 

[80] T. Shiromizu, Inverse cascade of primordial magnetic field in mhd turbulence, Phys. Lett. B 443 (1998) 127-130. 

[81] P. Schaeffer, On the dissipation rate coefficient in homogeneous isotropic decaying and forced turbulence, arXiv:1108.0575 

(2011) . 

[82] I. Proudman and W. H. Reid, On the decay of a normally distributed and homogeneous turbulent velocity field, Philosophical 
Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences (1954) 163-189. 

[83] G.L. Eyink and D.J. Thomson, Free decay of turbulence and breakdown of self- similarity, Physics of Fluids 12 (2000) 477. 

[84] P.D. Ditlevsen, M.H. Jensen and P. Olesen, Scaling and the prediction of energy spectra in decaying hydrodynamic turbulence, 
Physica A 342 (2004) 471-478. 

[85] G. Comte-Bellot and S. Corrsin, Simple Eulerian time correlation of full- and narrow-band velocity signals in grid-generated, 

"isotropic" turbulence, J. Fluid Mech. 48 (1971) 273 - 337. 
[86] Q. Chen, S. Chen and G.L. Eyink, The joint cascade of energy and helicity in three-dimensional turbulence, Phys. Fluids 15 

(2003) 361. 

[87] D. Bernard and K. Gawcdzki, Scaling and exotic regimes in decaying Burgers turbulence, J. Phys. A 31 (1998) 8735. 
[88] A. A. Fedorenko, P. Le Doussal and K.J. Wiese, in preparation. 

[89] P. Le Doussal and K.J. Wiese, 2-loop functional renormalization for elastic manifolds pinned by disorder in N dimensions, 

Phys. Rev. E 72 (2005) 035101 (R), cond-mat/0501315. 
[90] G. Blatter, M.V. Feigel'man, V.B. Geshkenbein, A.I. Larkin and V.M. Vinokur, Vortices in high-temperature superconductors, 

Rev. Mod. Phys. 66 (1994) 1125. 



Functional renormalization- group approach to decaying turbulence 



38 



[91] T. Giamarchi and P. Lc Doussal, Elastic theory of flux lattices in the presence of weak disorder, Phys. Rev. B 52 (1995) 1242—70. 
[92] L.D. Landau and E.M. Lifshitz, Fluid Mechanics, Volume 6, Pergamon Press, Oxford, 2 edition, 1987. 
[93] R.H. Kraichnan, Inertial-range transfer in 2- dimensional and 3- dimensional turbulence, J. Fluid Mech. 47 (1971) 525. 
[94] G. Falkovich and V. Lebedev, Universal direct cascade in 2-dimensional turbulence, Phys. Rev. E 50 (1994) 3883-3899. 
[95] J. Paret and P. Tabeling, Intermittency in the two-dimensional inverse cascade of energy: Experimental observations, Phys. 
Fluids 10 (1998) 3126-3136. 

[96] J. Paret, M.C. Jullien and P. Tabeling, Vorticity statistics in the two-dimensional enstrophy cascade, Phys. Rev. Lett. 83 (1999) 
3418-3421. 

[97] C. Pasqucro and G. Falkovich, Stationary spectrum of vorticity cascade in two-dimensional turbulence, Phys. Rev. E 65 (2002) 
056305. 

[98] D. Bernard, Three-point velocity correlation functions in two-dimensional forced turbulence, Phys. Rev. E 60 (1999) 6184-6187, 
chao-dyn/9902010. 

[99] D. Bernard, Influence of friction on the direct cascade of the 2d forced turbulence, Europhys. Lett. 50 (2000) 333, chao- 
dyn/9904034. 

[100] G.L. Eyink, Dissipation in turbulent solutions of 2D Euler equations, Nonlinearity 14 (2001) 787-802. 

[101] P. Constantin, On the Euler equation of incompressible fluids, Bulletin (new series) of the American Mathematical Society 44 
(2007) 603-621. 

[102] P. Constantin and F. Ramos, Inviscid limit for damped and driven incompressible Navier-Stokes equations in R 2 , 

Communications in Mathematical Physics 275 (2007) 529-551, arXiv:math/0611782. 
[103] P. Le Doussal and K.J. Wiese, Derivation of the functional renormalization group fS-function at order l/N for manifolds pinned 

by disorder, Nucl. Phys. B 701 (2004) 409-480, cond- mat/0406297. 
[104] P. Le Doussal and K.J. Wiese, Functional renormalization group at large N for disordered elastic systems, and relation to replica 

symmetry breaking, Phys. Rev. B 68 (2003) 174202, cond-mat/0305634. 
[105] P. Le Doussal and K.J. Wiese, Functional renormalization group at large N for random manifolds, Phys. Rev. Lett. 89 (2002) 

125702, cond-mat/0109204vl. 

[106] G. Boffetta, F. De Lillo and S. Musacchio, Inverse cascade in Charney-Hasegawa-Mima turbulence, Europhys. Lett. 59 (2002) 
687-693, nlin/0203044. 

[107] G. Falkovich and S. Musacchio, Conformal invariance in inverse turbulent cascades, arXiv: 1012.3868 (2012). 
[108] F. David, B. Duplantier and E. Guitter, Renormalization theory for interacting crumpled manifolds, Nucl. Phys. B 394 (1993) 
555-664. 

Contents 



1 Introduction 2 

2 Model and notations 3 

3 Known results and phenomenology 4 

3.1 Decaying Burgers 4 

3.2 Decaying Navier Stokes 5 

4 FRG equations 6 

4.1 Loop expansion: General strategy 6 

4.2 FRG equation for Burgers 7 

4.3 FRG equation for Navier-Stokes in momentum space 8 

4.4 FRG equation for Navier-Stokes in real space: isotropic turbulence 9 

4.5 Energy conservation and energy anomaly 9 

4.6 Enstrophy conservation and enstrophy anomaly in dimension 2 10 

5 Short-distance analysis, cusp or no cusp? 11 

5.1 Burgers 11 

5.2 Navier-Stokes 12 

6 Two-dimensional decaying turbulence (N = 2) 12 

6.1 Basic properties 12 

6.2 Isotropic turbulence 13 

6.2.1 FRG equations 13 

6.2.2 Searching for an isotropic fixed point 14 

6.2.3 Small-fc expansion 14 

6.2.4 Large-fc expansion 15 

6.2.5 Numerical solution 15 

6.2.6 Physical interpretation of the solution 16 

6.3 Periodic 2D-turbulcncc 16 



Functional renormalization- group approach to decaying turbulence 39 

7 Analysis of the FRG equation in three dimensions (N = 3), and in large dimensions (N — > oo) 18 

8 Decaying surface quasi-geostrophic turbulence 19 

9 Conclusions and Perspectives 20 
A 1-loop FRG equation in real space 20 
B 1-loop FRG equation in Fourier space 21 
C FRG equation for N = 3 periodic flows 23 
D Generating functional approach and diagrammatics 23 

E Distance geometry for FRG equation 26 

E.l Navier-Stokes 26 

E. 2 Burgers 27 

F Short-distance expansion of A(it): Amplitudes B and C for Burgers, Navier Stokes and surface 

quasi-geostrophic turbulence 27 

F. l Burgers 27 

F.2 Navier-Stokes 28 

F.3 Surface quasi-geostrophic turbulence 29 

G Small-fc expansion for Navier-Stokes in dimension N = 2 31 

H Asymptotic large-fc behavior for Navier-Stokes in N = 2 31 

I Numerical solution for the fixed point in dimension N = 2 33 

J Navier-Stokes equation in the limit of large N 34 



